^) 7 <£- //a// 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


Technical Memorandum 33-512 

Program Listing for Fault Tree Analysis 
of JPL Technical Report 32-1542 

Paul O. Che/son 



JET PROPULSION LABORATORY 

CALIFORNIA INSTITUTE OF TECHNOLOGY 
PASADENA, CALIFORNIA 


December 1 , 1 971 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


Technical Memorandum 33-512 

Program Listing for Fault Tree Analysis 
of JPL Technical Report 32-1542 

Paul 0 . Che/son 


JET PROPULSION LABORATORY 

CALIFORNIA INSTITUTE OF TECHNOLOGY 
PASADENA, CALIFORNIA 


December 1, 1971 



PREFACE 


The work described in this report was performed by the Quality 
Assurance and Reliability Division of the Jet Propulsion Laboratory. 


JPL Technical Memorandum 33-512 



CONTENTS 


I. Introduction 1 

II. Program Listing Z 

A. Main program FMAIN Z 

B. Subroutine FLT 13 

C. Subroutine READDS Z1 

D. Subroutine MODIFY ZZ 

E. Subroutine PRTEQ Z4 

F. Subroutine PRRBD Z5 

G« Subroutine IBSTBY Z6 

H. Subroutine SUB7 Z7 

I. Subroutine GATE Z8 

J. Subroutine SETIB Z9 

References 30 


JPL Technical Memorandum 33-51Z 


v 



ABSTRACT 


This technical memorandum presents the computer program listing 
for the fault tree analysis of Technical Report 32-1542, Reliability Compu- 
tation Using Fault Tree Analysis , Jet Propulsion Laboratory, Pasadena, 
Calif. , Dec. 1, 1971. The program is written in FORTRAN V and is cur- 
rently running on a UNIVAC 1108. 
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I. INTRODUCTION 


This technical memorandum presents the computer program listing 
for the MAIN program and those subroutines unique to the fault tree analy- 
sis described in Ref. 1. These subroutines are FLT, READDS, MODIFY, 
PRTEQ, PRRBD, IB STB Y, SUB7, GATE, and SETIB. The other sub- 
routines called by MAIN are used for analyzing the reliability block diagram 
of Ref. 1 and are listed in Ref. 2. The program is written in FORTRAN V 
and is currently running on a UNIVAC 1108. 
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TR£A*FTRE£.FMAIN 


1 

r 

FAULT TP5E COMPUTATION PROGRAM 

FM N 

10 

2 

c 


F MN 

20 

3 

c 


FM N 

30 

4 

c 

WRITTEN BY PAUL CHELSON 

F MN 

40 

5 

c 

JET PROPULSION LAb 

FM N 

50 

6 

c 

SECTION 153 — RELIABILITY 

F MN 

60 

7 

c 


FM N 

70 

3 

c 

LATEST EDITION - - 12 JULY 1971 

F MN 

80 

3 

c 


FM N 

9 C 

10 

c 


F MN 

ICC 

11 

c 

UN IV AC ilEfc* FORTRAN V VERSION* 65K CORE STORAGE REQUIRED. 

FMN 

lie 

12 

c 


F MN 

120 

13 

c 


•FM N 

130 

14 

c 


F MN 

140 

15 

c 

DATA DECK SETUP. 

FM N 

150 

lb 

c 


F MN 

160 

17 

c 

FAULT TREE DESCRIPTION. 

FMN 

17C 

18 

c 

IF THERE ARE N SASIC FAULTS* THERE WILL BE N CARDS OF THIS TYPE. 

F MN 

180 

19 

c 

BLOCK NUMBER (1-2) tFAULT PATH*3-621 WITH FORMAT 20A3* A 1 IN (80 

FMN 

190 

20 

c 

IN LAST CARD OF THIS TYPE. 

FMN 

200 

21 

c 


FMN 

210 

22 

c 

DISTRIBUTION TYPE (1-2)* MISSION TIMEC3-14) E FORMAT. - 

FMN 

22C 

23 

c 

1 - EXPONENTIAL* 2 = NOT AVAILABLE. 

FMN 

230 

24 

c 


FMN 

240 

25 

c 

ACTIVE PARAMETERS FCR EACH BLOCK OTHER THAN SENSE BL 0 CKS < SWI T CHES ) FM N 

250 

25 

c 

RG IS SET EQUAL TO I.C IF(15-25> LEFT BLANK. 

FMN 

26C 

27 

c 

ELCCK NUMBER (1-2 ) * LAMBDA (3-14 ) E FORMAT* R0C15-25) F FORMAT. 

FMN 

270 

23 

c 

LAST CARD IN THIS .SERIES HAS L AS T (80 SET EGUAL TO 5. 

FMN 

280 

23 

c 


FMN 

29 0 

30 

c 

IF NC DORMANCY INVOLVED. A BLANK CARD* OTHERWISE - 

FMN 

3CC 

31 

c 

&LANMI-2). DORMANCY FACTGR(S-14) E FORMAT. TO BE MULTIPLIED BY 

FMN 

310 

32 

c 

ACTIVE LAMBDAS OR S9(l-2)» 3 lANK<3-801 AND READ INDIVIDUAL OORMANTF MN 

320 

33 

c 

LAMEDAS - BLOCK NUMB ER ( 1-2 ) . LAMBDA DORMANT ( 3-14) E FORMAT. 

FMN 

33 G 

34 

c 

LAST CARD IN THIS SERIES HAS LAST (30 SET EQUAL TO 6. 

FMN 

34C 

35 

c 


FMN 

35C 

3 5 

c 


FMN 

36C 

37 

c 

SWITCHING CP II C NS FOR EACH SENSE BLOCK. 

FMN 

37 C 

‘ 38 

c 


FMN 

38G 

39 

c 

0 - NO SWITCHING l PROBABILITY OF SWITCH WORKING = 1.0). 

FMN 

390 

40 

c 

SENSE 3LCCK(l-2). 

FMN 

4 00 

41 

c 


FMN 

41C 

42 

c 

1 = CONSTANT PROBABILITY THAT SWITCH WORKS. 

FMN 

420 

43 

c 

SENSE BLOCK ( 1-2 > t&LANKt 3-14 ) .PR OE AblL IT Y ( 15- 2 5 ) F FORMAT. 1(80). 

FMN 

430 

44 

c 


FMN 

4 40 

45 

c 

2 = DORMANT FAIlURE RATE FOR SWITCH. 

FMN 

450 

46 

c 

SENSE 3L0CKC1-2) .LAMBDA DO R MA NT ( 3- 14) E FORMAT* 2(30. 

FMN 

4 60 

47 

c 


FMN 

47 C 

43 

c 

3 r DORMANT AND ACTIVE FAILURE RATE FOR SWITCH! 2 CARDS/SWITCH). 

FMN 

4 80 

49 

c 

SENSE BLOCK (1-2 ) * LAMBDA D ORM AN T ( 3 -14 ) l FORMAT* 3 (eC ) . 

FMN 

4SC 

50 

c 

SENSE BlCCK(1-2) .LAMBDA AC T IVE ( 3-14 I . 

FMN 

5 00 

51 

c 


FMN 

51C 

52 

c 

4 = NOT AVAILABLE. 

FMN 

5 20 

53 

c 

5 r NOT AVAILABLE. 

FMN 

530 

54 

c 


FMN 

5 40 

55 

c 

LAST CARD IN DATA DECK HAS LAST (60) SET EGUAL TO 7. 8. OR 9 - 

FMN 

550 

56 

c 

7 - RECALCULATE WITH NEW PARAMETERS* 3 z NEW DIAGRAM* 9 = END. 

FMN 

560 

57 

c 


FMN 

57C 

58 

C*****************************«**«**************************************fm\ 

5 80 

59 

c 


FMN 

590 

60 


COMMON /ALLSUB/M1*M2*M3*M4*M5* M6»M7*Li *L3»L4*L5*L& 

FMN 

SCO 

c 1 


CCMMGN/LGSUB/I6 (15*2 *50 ) * IIS (5G >* IR6 (5C* 3C ) • ITEMP(41 *2C ) »N1* 

FMN 

610 

62 


1 I S AVE ( 50.2GC ) .ISUC* JSUC *JSl* JSAVE (50* 200 • IJS1. T OT * PR OB ( 6 5 ) . 

FMN 

6 20 

63 


2 PR0BK&5 )*IFRINT .NSTDb Y( 15) »N SBYM X tNOPR NT 

FMN 

630 

64 


COM MGN / SM SAJ.3Z.fti R_UP-(.2.0C.L»N.EJE.MX-».N3 0U-T-(-5 0 )-.-N3-0T MX-.-N B-I-N-(-5G‘)'>-N£ rNMX^t 

"FMN" 

-640 

65 


1 N6NUM ( 50 ) • NbNM AX 

FMN 

650 

66 


DIMENSION RG15C>»NP(5C)fRT(65*lCi) * PSNS 1 1 D 1) * IS t 5 u ) * T L l 50 ) 

FMN 

6 60 

67 


DIMENSION IN (15 ) *IGUT (15). NRt ( 3C ) .KR B (30 *KS AVE ( 50 *2 DC ) 

FMN 

67 C 

68 


DIMENSION LSAVE(5C# 2CGI* TlD(53) 

FMN 

6 80 

69 


LIME NSICN SWPR0B(5L ) .TLCS (5G ). TLS( 50) 

FMN 

69 0 

7 C 


DOUBLE PRECISION DPR EC .DPR ECi * DPR E C2 . DPREC 3 . OPRE C4 . DPREC5 

FMN 

7 CO 

71 


DCUELE PR ECI SI GN PIN T * Y AV l tF SN S . TO T . PRO B .PRO B1 

FMN 

71C 
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72 

C 


FMN 

7 20 

73 

c 


FMN 

73 0 

74 

c 


FMN 

74G 

75 


Ml = 50 aKAXIMJM NUMBER OF BLOCKS IN DIAGRAM. 

FMN 

75 C 

7 6 


M2 = 200 E M AX I MUM NUMBER CF SUCCESS PATHS. 

FMN 

7 60 

77 


M3 = 14 EM A XI M UK NUMEE n CF INPUT S/G UTP U TS T 0/ F ROM GNE BLOCK. 

FMN 

77 C 

73 


M4 = 15 EMAXIMUM NUMBER OF SENSE BLOCKS CONTROLLING STANDS*. 

FMN 

7 8C 

7S 


ME r 29 DM A XI M UK NUMBER CF STANDBY BLOCKS CONTROLLED 6Y 1 SENSOR 

.FMN 

79 C 

30 


Mo = 20 i MAX IMUM NUMBER OF EQUIVALENT BLOCKS IN A SINGLE SET. 

FMN 

8CC 

31 


M 7 = 20 EH A XI MUM NUMBER GF EQUIVALENT BLOCK SETS. 

FMN 

810 

82 


Li=Ml*l 

FMN 

3 20 

63 


L3=M3*1 

FMN 

830 

64 


L4 - M4 ♦ Ml 

FMN 

84C 

£5 


L 5 = MS 

FMN 

65C 

86 


LS= ( 2 * Mo ) ♦ i 

FMN 

8 60 

67 


LL6-K6+1 

FMN 

87 0 

38 

c 


FMN 

8 83 

89 

c 


FMN 

69 C 

90 

c 

THE FOLLOWING SHOWS THE RELATIONSHIP OF THE ABOVE TO STORAGE. 

FMN 

9GC 

91 

c 


FMN 

SIC 

92 

c 

COMMON/ AlL SUB/ M i * M 2 # M 3 *M 4* H 5* Mb * M7 . Li * L3 * L 4 *L 5 * L6 

FMN 

920 

S3 

c 

CGMMCN/LGSUb/IbC L3.2 *H1 ) t IIS (Ml ) » 1 RE ( Ml t L5 J • I TEMP ( L6 t M7 ) »N1 • 

FMN 

930 

94 

c 

1 IS AVE( Ml » M2 ItlSUCt JSUC *JSlt JSAVEtMi* M2 > * I J3 It T OT * PR 03 ( L 4 ) * 

FMN 

940 

95 

c 

2 PRC5l(L4)*lFRINTtNSTDBY(M4) .NS6YMX 

FMN 

950 

96 

c 

COM MCN/ SM S U 3/ NP UP ( M2) »N PUP MX . No CUT ( Ml) t \S CTMX i NBINI Ml ) • NEINMX • 

FMN 

960 

97 

c 

1 KBNUKC Ml ) tNBNKAX 

FMN 

97 D 

98 

c 


FMN 

983 

S9 

c 

NS T - NSTEF ♦ 1 

FMN 

9SC 

IOC 

c 

DIMENSION RO(Ml)tNF (Kl)t RT(L4.NST) » PENS ( NS T ) i IS ( Ml ) » T L ( Ml ) 

FMN 

1COO 

1C1 

c 

DIMENSION IN (Ml ) tZCUTCKl) • NRE < L5 ) tKR s ( L5 ) tKS AVE t Ml. M2 ) 

FMN 

1C1C 

1C2 

c 

DIMENSION LSAVEtMlt M2 )• IDR IK1) • TLO (Ml) 

FMN 

1020 

103 

c 

0IMEN5I0N SWPRC b ( Ml ) • TL DS ( Ml ) t TL S ( Ml ) 

FMN 

103 C 

104 

c 


FMN 

1G 40 

1C 5 

c 


FMN 

1050 

ICS 

c 

CLEAR VARIABLES TO ZERO AND SET CONSTANTS. 

FMN 

1060 

1C-7 

c 


FMN 

1C 70 

1 C 3 

c 

IPRINT CONTROLS HCW MUCH OUTPUT WILL BE PRINTED. 

FMN 

1080 

109 

c 

IPRINT-C - PRINTS F L T TREE AND RESULTS. 

FMN 

10SC 

110 

c 

IPRINT - 1 - C ♦ PRINTS THE EQUIVALENT 3LCCK DIAGRAM. 

FMN 

nee 

111 

c 

IPRINT -2 - 1 + OVERALL SYSTEM PROBABILITY TREES ARE PRINTED. 

FMN 

me 

112 

c 

IPRINT- 3 - 2 + ALL PROBABILITY TREES ARE PRINTED. 

FMN 

112C 

113 

c 

IPRINT-4 - 3 ♦ ALL DIAGNOSTIC INF 0 IS PRINTED. 

FMN 

1130 

114 

c 


FMN 

1140 

115 


I F R I N T = 0 

FMN 

use 

115 

c 


FMN 

116G 

117 

c 

THE VARIABLE NSIG SPECIFIES THE NUMBER OF *S IGNIF IC ANT * FIGURES 

FMN 

1170 

118 

c 

TC 6E PRINTED FOR THE RELIABILITIES. NOTE THAT *5 1 GNIFICA NT* 

FMN 

1180 

119 

c 

FIGURES IS DEFINED AS THE NUMBER OF NON-NINES IN THE RELIABILITY 

FMN 

1 19 C 

120 

c 

NUMBER. THUSf .93986* .385* AND. 56 ALL HAVE TWO •SIGNIFICANT* 

FMN 

12CG 

121 

c 

DIGITS . 

FMN 

1210 

122 

c 


FMN 

1220 

123 


NS IG -3 

FMN 

1230 

124 

c 


FMN 

1240 

125 

_ c 

'THE ARR'A'Y ~V~ HOC DO T HE “F CR MAT _ F 0 R ~T HE - FTN A L: ‘R E'LrAEflLTT Y~ P RINT ING - ^" 

FMN" 

1250 - 

126 

c 

THE ARRAY VK HOLDS THE ALPHA NUMBERS NEEDED TO CHANGE THE 

FMN 

1260 

127 

c 

VARIABLE FORMAT V. 

FMN 

1270 

128 

c 


FMN 

1280 

129 


DIMENSION V ( ID ) »VK(1C) 

FMN 

1290 

130 

c 


FMN 

1300 

131 


DATA V / * f DbHLR* **EL1ABI***LI7Y 0*t*F THE *.*SYSTEM*** THRU *. 

FMN 

1310 

132 


l’TlME’ * •FIC.C * f * * 9H HG U* * * R S - * . • F • • * 1C • # * . * * * 3 * # f ) * / 

FMN 

1320 

133 


0 AT A VK/* 1* . *2 * . ’3' * ’ 4* .♦ 5* • *6 * t »7* * »8* • * 9' * *10 V 

FMN 

1330 

134 

c 


FMN 

1340 

135 

595 

L AST = 0 

FMN 

1350 

136 


NCPRNT-1 

FMN 

1360 

137 


DC 19 I- 1 »L3 

FMN 

1370 

138 


DO 19 J- 1 * 2 

FMN 

1380 

1 33 


DO 15 K = i .Ml 

FMN 

1390 

140 

19 

I5( I. J.K ) = Q 

FMN 

14 CC 

141 


DC 6 1=1. Ml 

FMN 

1410 

142 


IS( I)=Q 

FMN 

142C 

143 


II S ( I ) =D 

FMN 

1430 

144 


00 6 K = 1 * L 5 

FMN 

1440 

145 

6 

IRBCItKirC 

FMN 

145C 
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145 

C 




FMN 

146C 

147 


997 REWIND 1C 



FMN 

14 7 G 

143 


NSB YMX=G 



FMN 

1480 

149 


KT-C 



FMN 

149 G 

15C 


N6H0LD=D 



FMN 

1 5 CD 

151 


N5TEP-1GC 



FMN 

1510 

152 


JJ J=C 



FMN 

1520 

153 


KR 61 =1 



FMN 

153 C 

154 


NR31=G 



FMN 

154Q 

155 


DO 1C 1 = 1. L4 



FMN 

155 C 

156 


PRCBI I ) = C. ZQ+G 



FMN 

1560 

157 


PR061CI ) =C* C D + C 



FMN 

157 D 

isa 


00 1C L= L # 101 



FMN 

158C 

159 


RTII t L ) = C • 



FMN 

153 C 

16C 


10 PSNS(L)=G.GD+C 



FMN 

16C0 

161 


DC 15 1 = 1 » M 4 



FMN 

161 C 

162 


15 NST D3 Y ( I ) = G 



FMN 

16 20 

163 


DO 16 1 = 1. L3 



FMN 

163 C 

164 


INC 11=0 



FMN 

1640 

165 


16 IOUT II >=G 



FMN 

1650 

166 


DO 2C 1=1. Ml 



FMN 

1660 

167 


S WPR C6(2)=1»C 



FMN 

167 C 

168 


tldsi I) = C. 



FMN 

1680 

169 


TLS C I > =0 • 



FMN 

169 C 

17C 


tli I)=0. 



FMN 

17 QC 

171 


TLDI IJ=C. 



FMN 

1710 

172 


ROC Il=C. 



FMN 

1720 

173 


NF (I 1=0 



FMN 

1730 

174 


DO 17 K = 1 * L 5 



FMN 

1740 

175 


NR 6 1 K ) =0 



FMN 

1 75 C 

176 


17 K Rd ( K ) =0 



FMN 

176C 

177 


DC 20 K=1*M2 



FMN 

1770 

178 


JSAVEI I * K ) =C 



FMN 

1780 

179 


20 I S AVE (I * K ) = 0 



FMN 

1790 

18Q 

C 




FMN 

18C0 

161 


IF (LAST. £0.7 1 60 TO 3G0C 



FMN 

1610 

132 

C 




FMN 

1820 

183 

C 

READ FAULT TREE AND CONVERT TO BLOCK DIAGRAM. 



FMN 

lezc 

134 

c 




FMN 

184C 

165 


CALL FLTINE.NSRi 



FMN 

1850 

136 


RG C Ml 1 = 1.0 



FMN 

1860 

167 


TL(Kl) = C.L 



FMN 

1670 

188 


TLDIMl) = G.G 



FMN 

188C 

189 

c 




FMN 

iesc 

19C 

c 

N3H0LD = TEMPORARY HOLD FOR THE HIGHEST NUMBER 3 LOC Kt N8MAX ) . 

FMN 

1900 

191 


NBHOLD = NE 



FMN 

1910 

192 

c 




FMN 

1920 

193 

c 

ISUC = SUCCESS BLOCK NUMBERIEASE OF PROBABILITY TREE TO BE GENER- 

FMN 

153 C 

194 

c 

ATED BY CALLING SUBROUTINE TREE). 



FMN 

1940 

195 

c 




FMN 

1950 

196 


017 ISU C=N8 



FMN 

I960 

197 

c 




FMN 

197 C 

198 

c 

JSUC AND IPRINT ARE USED IN SUBROUTINE TREE TO CONTROL PR 

INTING. 

FMN 

19 80 

159 

c 




FMN 

199 C 

2CQ 


JSU C= NB 



FMN 

2GQ0 

201 

c 




FMN 

2010 

2C2 

c 

NBMAX = THE HIGHEST BLOCK NUM3ER CF SLOCKS 1- 5C/ 1- 

Ml. 


FMN 

2C20 

2C3 

c 




FMN 

203 C 

2C4 


NBM AX = N8HCLG 



FMN 

2040 

205 

c 




FMN 

2C5C 

2 C 6 

c 

SET THE FIRST ELEMENTS OF THE 13 ARRAY WITH THE QUANTITY 

OF 

FMN 

2C6Q 

2C7 

c 

INPUTS AND OUTPUTS TO EACH SlOCK. 



FMN 

207 C 

20 8 

c 




FMN 

2080 

2C9 


CALL SETIBdbtMl »L3 ) 



FMN 

2090 

210 — 

C“ 




FMN 

21C0 

211 

c 

READ SENSE BLOCK AND STANDBY BLOCKS OF THAT SENSE 

BLOCK. 


FMN 

2110 

212 

c 




FMN 

2120 

213 

c 

IIS ARRAY HOLDS SENSE BLOCKS OF CRIGINAL TREE. 



FMN 

2130 

214 

c 

IS ARRAY HOLDS SENSE BLOCKS OF ORIGINAL TREE AND 

VARIES. 


FMN 

2140 

215 

c 

NNSR = NUMBER OF SENSE BLOCKS IN ORIGINAL TREE. 



FMN 

215 C 

216 

c 

NS R = NUMBER CF SENSE BLOCKS IN CRIGINAL TREE AND 

VARIES 

* 

FMN 

2160 

217 

c 

IRB ARRAY HOLDS STANDBY BLOCKS CONTROLLED BY THEIR 

SENSE 

BLOCKS. 

FMN 

2170 

218 

c 

I RB ( NUMBER OF SENSE BLOCK. 1) = NUMBER OF 

STANDBY 

BLOCKS 

FMN 

218C 

219 

c 

CF THAT SENSE 6 LOCK • 



FMN 

21SC 

220 

c 

I RB ( NUM3ER OF SENSE BLOCK. J > = STANDBY BLOCK NUMBER CONTROLLED 

FMN 

2200 

221 

c 

BY THAT SENSE BLOCK* J = 2 • • • * 



FMN 

221 G 
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222 

C 


FMN 

22 20 

223 


NN SR -NSR 

FMN 

223 C 

224 


IF(NSR.£Q.GI GO TO 84Q 

FMN 

2240 

225 


DO 2C01 1=1. NSR 

FMN 

225 G 

228 


IS( I> = IISC ) 

FMN 

2260 

227 


NS = IIS(1 > 

FMN 

2270 

228 


00 2CC5 J = L 5 . 2 . — 1 

FMN 

228C 

229 


I R 6 ( NS . J ) = IRB(NS.J~1> 

FMN 

2290 

230 

2 COS 

CONTINUE 

FMN 

23C0 

231 


DO 2GC2 J = 1 » M5 

FMN 

2310 

232 


IF IIR3(N3.J+1> .EQ.d GO TO 2003 

FMN 

2320 

233 

2 GO 2 

CONTINUE 

FMN 

2330 

234 

2003 

I R3 ( NS • 1 ) - J-l 

FMN 

2340 

235 

2CG1 

CONTINUE 

FMN 

2350 

2 3 S 

C 


FMN 

2360 

237 

C 

standby llccks must follow all OTHER INPUT BLOCKS IN INPUT LISTS 

FMN 

237C 

238 

C 

FOR 12. ARRANGE 19 AND PRINT THE SENSE 3L0CKS AND THE STANDBY 

FMN 

2380 

2 3S 

C 

ELOCKS. 

FMN 

2390 

2 4 C 

C 


FMN 

24 Cu 

241 


CALL IBST6Y { JS28 »NSR .IS rl RE r IB . IPRIN 7 » L3 t L5 ) 

FMN 

2410 

242 

c 


FMN 

2420 

243 

c 

READ AND GENERATE ITEMF ARRAY OF EQUIVALENT ELOCKS. 

FMN 

243 G 

244 

c 

N1 = MAXIMUM NUMBER OF EQUIVALENT BLOCK SETS. 

FMN 

2440 

245 

c 


FMN 

245 C 

246 

c 


FMN 

24 60 

247 

c 

PRINT RELIABILITY BLOCK DIAGRAM. 

FMN 

247 C 

248 

c 


FMN 

24 8C 

249 

640 

IF (IPRINT.EG.C ) 60 TO 2C45 

FMN 

2490 

25D 


CALL PRRSD ( 13 » Ml ) 

FMN 

2500 

251 

c 


FMN 

251C 

252 

c 

PRINT EQUIVALENT BLOCKS. 

FMN 

2520 

253 

c 


FMN 

253 C 

254 

2045 

IF (Nl.EQ.C) GO TO 2040 

FMN 

2540 

255 


CALL PR TEG (I TEMFtIFRiNT .N1.LL6 1 

FMN 

255 C 

25b 

2040 

CONTINUE 

FMN 

2560 

257 

C 


FMN 

257 C 

250 

c 

RELIABILITY BLOCK DIAGRAM INPUT CHECKED. 

FMN 

2580 

259 

c 


FMN 

259 C 

26C 

c 


FMN 

2600 

261 

c 

GENERATE ORIGINAL TREE. 

FMN 

2610 

262 

c 

ISAVE HOLDS THE PROBABILITY TREE SUCCESS PATHS 

FMN 

2620 

263 

c 

JS1 - NUMBER OF SUCCESS PATHS. 

FMN 

263 C 

264 

c 


FMN 

2640 

265 

3CCC 

CALL TR EC ( S 9 58 ) 

FMN 

265 C 

266 

c 


FMN 

2660 

267 

c 

DUPLICATE ORIGINAL TREE INTO JSAVE . 

FMN 

267C 

268 

c 


FMN 

2680 

269 


CALL DUPTRE< uSAVE.lJSitlSAVE » JS1) 

FMN 

2620 

27C 

c 


FMN 

2700 

271 

c 

ID IS T DETERMINES WHAT FAILURE DISTRIBUTION IS USE D ( 1 = E XP ONEN T I AL ) 

FMN 

2 71 C 

272 

c 

TTOT - MISSION TIME. 

FMN 

2720 

273 

c 


FMN 

2730 

274 


READ 40CG * IOIST. TTOT 

FMN 

2740 

275 

400 C- 

F-GRM A-UI 2-t El 2. 7J 

FMN 

2 75 0 

276 

c 


FMN 

27 60 

277 

c 

READ PARAMETERS OF FAILURE DISTRIBUTION. 

FMN 

2770 

278 

c 

RC IS SET = i.C IF READ IN AS BLANK OR ZERO. 

FMN 

2780 

279 

c 


FMN 

2720 

280 


GO TC <4100.4200. 4300.4400. 4SCG) .IDIST 

FMN 

23GC 

281 

c 


FMN 

2810 

282 

c 

READ ACTIVE PARAMETERS. 

FMN 

2820 

283 

c 


FMN 

2830 

284 

4 10C 

READ 4 1C 1 * Nc.TLAMbD. RC1.LA5T 

FMN 

2840 

285 

4 1C 1 

FORM AT ( 12 t El 2. 7 •FiC.7»5 5X .11 ) 

FMN 

285 0 

286 


IF< RG1141C9.411u.41C9 

FMN 

286C 

287 

4109 

TLCNB)=TLAMBD 

FMN 

2870 

288 


RC( N3 URQ1 

FMN 

2880 

289 


GO TC 412G 

FMN 

2890 

290 

4 1 1C 

RC( NS 1 = 1.0 

FMN 

2900 

291 


TL<NB>=TLAMED 

FMN 

2910 

292 

4120 

IF{ LAST.LT.5) GO TO 41CG 

FMN 

2920 

293 

C 


FMN 

2230 

294 


IFCNSR.EQ.O) GO TO 47 1 C 

FMN 

2940 

295 

c 


FMN 

2950 

296 

c 

READ DORMANCY PARAMETERS AND SWITCHING OPTIONS 

FMN 

2960 


JPL Technical Memorandum 33-512 


5 



2S7 

C 


FMN 

297 C 

298 


CALL REA003(N3Rf TLDt TL .R Of S WP R 03 # TLDS * TLS . IRE t IS > 

FMN 

2980 

2S9 

C 


FMN 

299 C 

300 

C 


FMN 

3000 

3C1 

C 

READ LAST(£C)» LAST = 7. 8* CR 9. 

FMN 

301 D 

3C2 

C 

7 = RECALCULATE WITH NEW PARAMETERS* 8 = NEW DIAGRAM. 9 - END. 

FMN 

302C 

303 

C 


FMN 

3D3C 

304 

4 7 1C 

READ 27C 7 * LAST 

FMN 

304C 

3C5 


GO TG 4800 

FMN 

2C50 

306 

C 


FMN 

3060 

307 

42C 0 

CONTINUE 

FMN 

307 C 

308 

4300 

CONTINUE 

FMN 

3G8G 

3G9 

4 4GC 

CONTINUE 

FMN 

3C9C 

310 

4500 

CONTINUE 

FMN 

3100 

311 

C 


FMN 

3110 

312 

c 

42GD.43CGf4400.45uC AVAILABLE IF YOU WANT TO USE OTHER THAN 

FMN 

3120 

313 

c 

EXPONENTIAL DISTRIBUTION. ADD APPROPRIATE STATKENTS FOR NEW 

FMN 

313 C 

314 

c 

DISTRIBUTION. 

FMN 

3140 

315 

c 


FMN 

315 0 

316 

4800 

IF(NSR.EQ.C) GO TO 8000 

FMN 

3160 

317 

c 


FMN 

3170 

318 

c 

FIND A SENSE TREE WITH NO OTHER SENSE BLOCKS IN IT. 

FMN 

3180 

319 

c 


FMN 

319C 

320 

700C 

CONTINUE 

FMN 

3200 

321 


DO 7010 JJ-1.K1 

FMN 

3210 

322 

7010 

NP ( JJ 1=0 

FMN 

3220 

323 


I NDE P = C 

FMN 

3230 

324 


M=1 

FMN 

3240 

325 


J=1 

FMN 

3250 

328 


1 = 1 

FMN 

32 60 

327 

7011 

IF ( JSAVE(I.J).NE.C) GO TO 7001 

FMN 

3270 

328 


IF ( 1-1 ) 7002. 70C3.70C2 

FMN 

3280 

329 

7GC1 

K = 1 

FMN 

329 C 

330 

7004 

L = 1 

FMN 

3 3 CO 

331 

700 5 

CONTINUE 

FMN 

3310 

332 


IF< NP(L ) .El.lS(K)) GO TO 7006 

FMN 

3320 

333 


IF (L.EQ.M) GOTO 7007 

FMN 

3330 

334 


L=L*1 

FMN 

3340 

335 


GO TO 7 CG5 

FMN 

335 C 

336 

7007 

XFCABSC JSAVEC I . J )) . El . IS (K ) > GO TO 7GC8 

FMN 

3360 

337 

7006 

IFCK.EQ.NSR) GO TO 7CC9 

FMN 

3370 

338 


K= K + l 

FMN 

3380 

339 


GC TO 7004 

FMN 

339 C 

340 

70C8 

IF( IS(KJ.EO.INDEF) GO TO 7003 

FMN 

3400 

341 


M=M*1 

FMN 

3410 

342 


NP< M— 1 ) = INDEP 

FMN 

3420 

343 


I NDE P = I S C K 1 

FMN 

343 G 

344 


IF ( IP RINT • NE • 4 ) GO TO 7CC3 

FMN 

3440 

345 


PRINT 2 DOC 1 . INDEP 

FMN 

3450 

346 

7009 

IF ( I.EQ.M1 ) GO TO 7002 

FMN 

3460 

347 


1=1*1 

FMN 

347C 

348 


GO TO 7011 

FMN 

34 80 

349 

7002 

IF ( J • E Q » M2) SO TO 7003 

FMN 

3490 

350 


1=1 

FMN 

3500 

351 


J= J*1 

FMN 

3510 

352 


GO TO 7011 

FMN 

3520 

353 

7003 

I SUC = 1 NDE F 

FMN 

3530 

354 


IF ( IP RINT. NE. 4) GO TO 7030 

FMN 

3540 

355 


PRINT 2CCC1. 1SUC 

FMN 

355 C 

356 

2CCC1 

FORMAT (13) 

FMN 

3560 

357 


PRINT 20CC2. (IS(I)tI=i .Ml) 

FMN 

357 C 

358 


PRINT 2000 2 . ( NP ( I ) * 1= 1. Ml ) 

FMN 

3580 

359 

2000 2 

FORM AT ( IHC25 (13 *1X)/1H 25 (13 .IX) ) 

FMN 

359C 

360 

C 


_fLMN_ 

-3600- 

381 

c 


FMN 

361 C 

362 

c 

DETERMINE WHAT TYPE OF STANDBY IS INVOLVED. 

FMN 

3620 

363 

c 


FMN 

3630 

364 

c 


FMN 

3640 

365 

c 

IF SENSE BLOCK CGNTRCLS ONE STANDBY BLOCK. SET ISTDBY = 1. SEE IF 

FMN 

365 G 

366 

c 

THAT STANDBY SLOCK HAS AN INPUT. IF NO INPUT GO TO 7700. IF IT 

FMN 

3660 

367 

c 

HAS AN INPUT GC TC 7500. IF SENSE BLOCK CONTROLS MOPE THAN ONE 

FMN 

367 C 

368 

c 

STANDBY BLOCK. SET IST0BY=2 AND CALL SUBROUTINE S T 0 BY 2 . 

FMN 

36 80 

369 

c 


FMN 

369 C 

370 

7030 

IF<IRti(lNO£P*i).GT.l) GO TO 7050 

FMN 

3700 

371 


I STD B Y = 1 

FMN 

3710 

372 


IRBB=IR3( INDEP. 2) 

FMN 

37 2C 
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373 


ITEST-IBdtl tlRBB ) 



FMN 

373C 

374 


IF< ITEST.EQ.C ) GO TO 77CC 



FMN 

3740 

375 


GO TC 7500 



FMN 

3750 

376 

7050 

1ST C3 Y-2 



FMN 

3760 

377 

C 




FMN 

377 G 

378 


CALL 3T0BY2(i3SSfINDEP) 



FMN 

3780 

379 

c 




FMN 

37SC 

380 

c 

DUPLICATE STANDBY TREE RETURNED IN ISAVE INTO LS AV E • 


FMN 

3 8 CO 

381 

c 




FMN 

3810 

302 


CALL DUPTREfL SAVE tLJSitl SAVE# JSl) 



FMN 

3320 

383 


IF (IPRINT.LE.2) 60 TO 7200 



FMN 

383 C 

384 


CALL TRPRNK2 » JSAVE.IJSl # 1 ) 



FMN 

3840 

385 


CALL TRPRNT (4#LSAVE»LJS1» 2 ) 



FMN 

365C 

386 

c 




FMN 

3860 

387 

c 

GENERATE FULL STANDEY TREE. 



FMN 

387C 

363 

c 




FMN 

38 80 

389 

c 

SUBROUTINE TRENUM STORES IN KBNUM ALL 

THE DIFFERENT 

ABSOLUTE 

VALUEFMN 

3890 

390 

c 

BLOCK NUMBERS IN THE STANDBY TREE. 

NUMOUT ZEROES OUT THESE 

FMN 

3900 

391 

c 

BLOCK NUMBERS. 



FMN 

3910 

392 

c 




FMN 

3920 

393 

7200 

CALL TRENUK(lSAVL#L JS i) 



FMN 

393C 

394 

C 




FMN 

3940 

395 


I S UC - 1 R 5 { I N D EP • 2 ) 



FMN 

3950 

396 

C 




FMN 

3960 

397 

C 

IF STANDBY BlOCK HAS NO INPUT# CC TO 

7500. 


FMN 

3970 

393 

C 




FMN 

3980 

399 


ITES7=IB<1#1 #ISUC ) 



FMN 

3990 

4 CO 


IF( ITEST.EQ.C ) GO TO 7 5C 0 



FMN 

4000 

401 

c 




FMN 

4010 

4 C 2 


CALL TR£E< $993) 



FMN 

4020 

403 

c 




FMN 

4C3C 

4C4 

c 

ZERO OUT STANDBY BLOCKS IN STANDBY TREE. 


FMN 

4Q4C 

4C5 

c 




FMN 

4050 

4C6 


CALL NUMOUT C$996#ISAV£» JSl) 



FMN 

4C60 

4C7 


CALL SFLEFT ($998 *ISAVEt USD 



FMN 

4D7C 

408 


CALL 0UPGUT(S993#ISAVE» JSl) 



FMN 

40 80 

4C9 


IF (IPRINT.LE.2) 60 TO 7400 



FMN 

4090 

410 


CALL TRPRNT (1 # IS AVE # JSl# 3 ) 



FMN 

41C0 

411 

7 4C0 

CALL TRENUHC1SAVE* JSl) 



FMN 

4110 

412 

C 




FMN 

4120 

413 

C 




FMN 

413 C 

414 

C 

CALCULATE THE PRO B ABIL IT Y ( F I TREE) THAT YOU NEED THE 

STANDBY 

TREE. FMN 

41 4G 

415 

c 




FMN 

415 G 

416 

c 




FMN 

4160 

417 

7 5CC 

I SUC - 1 NDE R 



FMN 

4170 

418 


CALL TREE ( $938) 



FMN 

4180 

419 

C 




FMN 

4190 

420 

C 

DUPLICATE SENSE TREE INTO KSAVE. 



FMN 

4200 

421 

c 




FMN 

4210 

422 


CALL DUPTR E (KSAVE»KJSl#ISAV£# JSl) 



FMN 

4220 

423 

c 




FMN 

4230 

424 

c 

SUBTRACT FULL STANDBY TREE FROM SENSE 

TREE GIVING PI 

TREE. 

FMN 

4240 

425 

c 




FMN 

4250 

A2 5. 


GO TO ( 7 515 #7 513 )_• ISTDBY 



FMN 

4260 

427 

c 




FMN 

427 C 

428 

c 

GENERATE FULL STANDBY TREE. 



FMN 

4280 

429 

c 




FMN 

429 C 

430 

7515 

ISUC=IR3( INDE P # 2 ) 



FMN 

4300 

431 


CALL TREE ( $9 5£ ) 



FMN 

4310 

432 

C 




FMN 

4320 

433 


CALL TRENUM { IS A VE * JSl) 



FMN 

433 C 

434 

C 




FMN 

4340 

435 

C 

INITIALIZE 



FMN 

435 C 

436 

C 




FMN 

436G 

437 

7518 

DG 7520 I - 1 » Ml 



FMN 

4370 

438 


NBOUT C I ) — G 



FMN 

43 80 

439 

7520 

CONTINUE 



FMN 

4390 

440 

C 




FMN 

4400 

441 

C 

ELIMINATE PATHS IN SENSE TREE IF THEY 

CONTAIN THE STANDBY BLOCKS. FMN 

441C 

442 

C 




FMN 

4420 

443 


NEGTKX-IREtINDEP #1 ) 



FMN 

4 43 C 

444 


DO 7530 J-i»N3CTMX 



FMN 

4440 

445 


NBOUT (J)-IRB(INDEP»J+1) 



FMN 

4450 

446 

7530 

CONTINUE 



FMN 

44 SO 

447 


CALL PATH0T($998. KSAVE. KJS1) 



FMN 

4470 

448 


IF< IPRINT.i_E.2> GO TO 75 35 



FMN 

4480 
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449 


CALL TRPRNT (3*KSAVE* KJS 1 1 4 ) 


FMN 

4490 

4 50 

C 



FMN 

4500 

451 

7535 

N6NMAX=NBNMAX+1 


FKN 

<*51 G 

452 


NBNUM ( NSNM AX) - IN DEP 


FMN 

4520 

453 


CALL NUMCUT($993»KSAVE»K JSl) 


FMN 

4530 

454 


CALL SF L EF T < $998*KSAVE»KJSl) 


FMN 

4540 

455 


CALL OUPOUT ($993*KSAV£*KJS1) 


FMN 

4550 

456 


IF (IPRINT.EG.G) GO TO 7538 


FMN 

4560 

457 


CALL TRPRNt(3»KSAV£fKJSl • 5 ) 


FMN 

4570 

458 

C 



FMN 

458 C 

459 

C 

CHECK PI TREE FOR STANDBY TREES REPLACEO 3Y STAND3Y BL OCK ( 51- 6 5 ) / 

FMN 

4590 

460 

C 

(L1-L4 l . 


FMN 

4600 

461 

C 



FMN 

461C 

462 

7538 

IF < NSBYMX.LE.l ) GC TO 7545 


FMN 

4620 

463 


REWIND 10 


FMN 

4630 

464 


DC 7 54 C N-I * NS6 YMX 


FMN 

4640 

465 


READl 10 ) NSStJSl 


FMN 

4650 

466 


DC 7537 I~1*M2 


FMN 

4660 

467 


READ(IO) ( ISAVE(J*I)» J-l * M 1 1 


FMN 

4670 

468 

7537 

CONTINUE 


FMN 

468C 

469 


CALL TRI,NTR($ 393 * K SAVE *K JS 1 #1 S AVE * JSl t NSTO BY* NS 9 YMX *NSS) 

FMN 

46 30 

4 70 

7540 

CONTINUE 


FMN 

47CC 

471 

C 



FMN 

4710 

472 

c 

DUPLICATE KSAVE INTO ISAVE. 


FMN 

4720 

473 

c 



FMN 

4730 

474 

7545 

CALL DUPTREdSAVE* JSlt KSAVE tKJSl ) 


FMN 

4 74 C 

475 

C 



FMN 

4750 

476 

C 

PI TREE NOW IN ISAVE. 


FMN 

4760 

477 

C 



FMN 

4770 

478 


IF CIPRINT.LE.2 ) GO TO 5CCC 


FMN 

4 78 C 

479 


CALL TRPRNTI1 tlSAVEt JS1 » 6 ) 


FMN 

4790 

480 


GO TO 5 COO 


FMN 

4800 

481 

C 



FMN 

4810 

482 

C 

GENERATE SENSE TREE. 


FMN 

482 C 

483 

c 



FMN 

4830 

484 

7700 

I SUC -I NCEP 


FMN 

4 84 C 

485 


CALL TREE! $998) 


FMN 

485C 

486 

C 



FMN 

4860 

487 

c 

INITIALIZE 


FMN 

4870 

468 

c 



FMN 

4 88 C 

439 


00 7710 I=1#M1 


FMN 

48 90 

4 SD 


NB OU T ( I ) -C 


FMN 

490 C 

491 


N3NUM d ) -0 


FMN 

4910 

4 92 

7710 

CONTINUE 


FMN 

492 C 

433 

C 



FMN 

4930 

4 94 

C 

ELIMINATE PATHS IN SENSE TREE IF THEY CONTAIN THE 

STANDBY BLOCK. 

FMN 

494 G 

495 

C 



FMN 

4950 

4 96 


NBCUT (1)-IRB(INDEP»2 ) 


FMN 

4960 

497 


NB0TMX=1 


FMN 

4970 

498 


CALL PATHGT ( $998 • IS A V E * JS1) 


FMN 

498 C 

439 


IF( IPRINT.LE.2) GO TO 77 20 


FMN 

4990 

5C0 


CALL TRFRNT ( 1»ISAVE » JSlt 7 ) 


FMN 

5CCC 

5C 1 

C 



FMN 

5010 

502 

C 

ZERO OUT SENSE BLOCK NUMBER IN SENSE TREE. 


FKN 

502 0 

503 

C 



FMN 

5030 

504 

7720 

NbNUKC 1>=INDEP 


FMN 

5C4C 

5C5 


N9NMAX=1 


FMN 

5050 

506 


CALL NUMOUT ( $998*1SAVE* JSD 


FMN 

5060 

507 


CALL SFLEFT ($393# ISAVE* JSl ) 


FMN 

5070 

508 


CALL DUPOUT ($998 • ISAVE* JSl) 


FMN 

5080 

509 


IF ( IPRINT .LE* 2) GO TO 5G 00 


FMN 

5090 

51C 


CALL TRPRNT ( 1 » I SAVE * JSlt 8 ) 


FMN 

51CC 

511 

C 



FMN 

5110 

512 

c 

CALCULATE PI* THE PROBABILITY YOU NEED SENSE TREE* 

STORE IN PSNS. 

FMN 

512 C 

5T3 

c 



~FMN~ 

51'30~ 

514 

5C0G 

T-0. 


FMN 

514C 

515 


NOW-O 


FMN 

5150 

516 


K-C 


FMN 

5180 

517 

9000 

CONTINUE 


FMN 

5170 

518 


GO TO (51 CC* 52 00*5300*5400*5 SCO )* IDIST 


FMN 

5180 

519 

C 



FMN 

5190 

520 

5100 

DO 5101 I -1 * NBM AX 


FMN 

5200 

521 


DPREC=— T *TL ( I ) 


FMN 

5210 

522 


PR0B1(I)-RG ( I ) * DEXP ( DPR EC ) 


FMN 

5220 

523 

5101 

CONTINUE 


FMN 

5230 

524 

C 



FMN 

524 0 
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KRB ARRAY CONTAINS STANDBY BLOCKS THAT HAVE HAD THEIR PROBABILITY FMN 
PREVIOUSLY CALCULATED. FHN 

- - - - - - - -- - FMN 

DC 5 1 1C L J = 1»KRS1 FMN 

IF( KR3( JJ.EQ.C) GO TO 5110 FMN 

NB-KRE(J) FMN 

PROBlt N3)-RT( N8tK+l ) FMN 

CONTINUE FMN 

FMN 

SUBROUTINE SYSP CALCULATES THE SYSTEM PROBABILITY - RETURNS TOT. FMN 

FMN 

CALL SYSP FMN 

FMN 

IFtNCW.E3.lJ GO TO 2001 FMN 

K-K + I FMN 

Tr <TTGT*K J/NSTEF FMN 

PSNSC K ) -i • 00+ 0 - TOT FMN 

FMN 

PREVENT DIVISION 3Y ZERO WHEN 13 SIGNIFICANT FIGURES LOST. FMN 

FMN 

IFt PSNSt KJ.GT .O.CD + O) GO TC 5160 FMN 

IF (K.EG.1J GO TO 5160 FMN 

FMN 

NEXT STATEMENT ELIMINATES DIVISION BY ZERC LATER WHEN CALCULATING FMN 
RT. THE NUMBER USED RESULTS FROM SUBTRACTING .999338999993999999 FMN 
FROM 1 .30 + 0 JUST 3EFCRE LOSING ALL SIGNIFICANT FIGURES. THIS WILLFMN 
CAUSE AN INTEGRATION ERROR IN CALCULATING RTf THE RELIABILITY OF FMN 
A SENSE BLOCK OR SENSE TREE. THE SYSTEM RELIABILITY MAY NOT EE FMN 
EFFECTED IF IT APPROACHES 1.0. FMN 


IF (K.EG.l J GO TO 5160 


PSNS(K) = +.1734723475 57 53C7C9D— C17 
IF (K.LE.NSTEF) GC TO 510C 
I ST £P = NST£P + i 

IF 1 1 PRINT * NE .4 ) GO TO £ CD C 
PRINT 5130# ( F SNS ( X ) » K-l t IS TEP ) 

FORM AT15D26.13 J 

GO TO 6000 

CONTINUE 

CONTINUE 

CONTINUE 

CONTINUE 

COMMENT STATEMENT AFTER 45CC. 

GO TO SOCG 


FMN 

FMN 

FMN 

FMN 

FMN 

FMN 

FMN 

FMN 

FHN 

FMN 

FMN 

FMN 

FMN 

FMN 

FMN 

FMN 


FMN 

NUMERICALLY INTEGRATE THE RELIABILITY OF THE STAND8Y BLOCKS. FMN 

TIME VARIABLE = JT = C-MISSICN TIMECTTOT). KT = 1-1G1 INCRE MENT S . FM N 

FMN 

FMN 

GO TO <6100#64Ca> #ISTD6Y FMN 

FMN 

METHOD FOR ONE STANDBY 0 LO C K( 1ST D3Y = 1 ) . FMN 

FMN 

NS=IR3< INDEP# 2) FMN 

RTCNB.lJrRO (NS J FMN 

-GC- -TO — ( S1G 2-#6 202 # 6 3 G2-* 64 C2 #-S5 02 ) »-ID 1ST F-MN- 

FMN 

ISTEP-NSTEP+i FMN 

DO 6170 KT=2tISTEP FMN 

DPREC = (-TTOT* (KT-1J/NSTEPJ *TLtNB) FMN 

DPREC1= (-TTGT* (KT-1)/N STEP l+TLS (INDEP) FMN 

PINT= R0(N3)*PSNS(i)*DEXP (OPR EC) *SWPR03( INDEP ) +0 EX P ( DPREC1 J FMN 

JTEND-K T-l FMN 

DC 6150 JT-1#JTEND FMN 

DPREC2r( (-TTGT*( KT-1 J/NSTEP) *TL (NB J) + ( (-TTOT* (JT-1J/NSTEP >* FMN 

1 (TLO(NS)-TL(NB))) FMN 

DPRE C3- 1 (-TTCT* (KT-1 l/N ST EP ) *TL (N6 ))+<(- T TOT ♦ ( JT J /NSTEP )♦ FMN 

L (TLD(N3)-TL( NB> > > FMN 

DPREC4-( (-TTOT* (KT-1 ) /NSTEP ) *TL S ( IND EP )>♦((- T T OT* ( JT -1 ) /NSTE P ) * FMN 

L ( TLD5( INDEP) ~TLS ( INDEP) )) FMN 

DPRE C5- ( ( -TTCT* (KT-1 J/NSTEP ) *TL S ( INDEP > ) + ( t- TTCT+ (JT J/NSTEP )• FMN 
L ( TLDS ( INDEP) -TlS ( INDEP ) )) FMN 

YAVE*( ( DEXP ( DPREC2 )*DEXP( DPREC4 ) ) + ( D EXP ( DPR E C3 ) * DEXP (DPREC5 ) ) J/2DCFM N 
PINT - PINT + R0(NB)*(PSNS(JT+1>-PSNS(JT>) ♦ YAVE * SW PROB ( INDEP ) FMN 
CONTINUE FMN 

RT( NBtKT) -PINT/P SNS (KT) FMN 

CONTINUE FMN 
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£01 


IF< IPRINT.NE.4) GO TO 60 10 

F MN 

6010 

6C2 


PRINT 5388* (RT ( Ns * K T ) * KTri.ISTEPJ 

FMN 

6C2C 

SC 3 

5338 

F0RMAT(8£i5.3) 

F MN 

60 30 

604 


GC TC 6010 

FMN 

604 C 

60 5 

6202 

CONTINUE 

FMN 

6050 

6C6 

6302 

CONTINUE 

FMN 

6C6C 

SC 7 

6402 

CONTINUE 

FMN 

6070 

SC 8 

6502 

CONTINUE 

FMN 

608 C 

6C3 

C SEE 

COMMENT STATEMENT AFTER 4500. 

FMN 

6G90 

SIC 

C 


FMN 

610 C 

611 

C 


FMN 

6110 

S12 

C 

METHOD FOR MORE THAN ONE STANDBY BLCCK(ISTCBY-2) 

FMN 

6 12 C 

613 

c 

RO FOR ALL STANDBY 3L OCK S < 5 1“ 65 > - 1 

FMN 

6130 

S14 

c 


FMN 

6140 

615 

6400 

N3-NSTD3Y( NS8YMX) 

FMN 

6150 

616 


R T (No*l ) -1.0 

FMN 

6160 

617 


GO TO (6103*6203*5303* 64 03* 65 03) *10 1ST 

FMN 

6170 

618 

c 


FMN 

6180 

619 

6103 

ISTEP=NSTEP+1 

FMN 

6190 

620 


DC 6 5C0 K T-2 iISTEP 

FMN 

6200 

621 


OPREC - ST3YPR(LSAVE*LJSlfRC*TL*RT»KT*0*TT0T* NST£P#TLO) 

FMN 

6210 

622 


DPRE Cl- (-7 TOT* t KT-D/NSTEP >*7LS t INOEP) 

FMN 

£220 

623 


PJNT = PSNS( 1 ) * DPREC* SWPRC 3( IND EP) ♦ DE XP ( DPREC 1) 

FMN 

6230 

624 


JTE ND-KT-1 

FMN 

6240 

625 


DO 6450 JT=1#JTENC 

FMN 

6250 

626 


0PREC4 = ( ( -TTGT* ( KT-1 )/NSTEP ) *TLS< XNOEP ) ) + ( ( -TTuT* (JT-1)/NSTEP ) * 

FMN 

626 C 

627 


1 (TLDS(INDEP)-TLS(INDEP) )) 

FMN 

6270 

628 


DPRE C5 = < < -TTGT*(KT-1)/NSTEP ) *TL S < INDEP ) ) ♦ C (-T TOT* (JT l/NSTEP )♦ 

FMN 

6280 

629 


1 < TLD5( INO£P)-TLS( INOEP) )) 

FMN 

6290 

63C 


J JT- JT-I 

FMN 

6 30 C 

631 


YAVE=l ST3YPR1 LSAVE.LJS1. RO.TL.RT.KT. JT . TT OT .NS T £ P. T LC J »OEX P < DPREC5F MN 

631C 

622 


1 ) ♦ 5TEYPR ILSAVE.LJSI . RO.TL.RT.KT. JJT.TTOT .NSTEP.TtDI *OEXP (CPRECFMN 

632 C 

633 


24 ) ) /2 m GO* 0 

FMN 

6 3 30 

634 


PINT r PINT •» (PSNSt J T+ 1 > -F‘ S NS < JT >) * YAVE * S WPR OB ( INOEP ) 

FMN 

6340 

635 

6450 

CONTINUE 

FMN 

6350 

636 


RT (Nfc*KT)-PIM/PSNS( KT) 

FMN 

6360 

537 

6500 

CONTINUE 

FMN 

6370 

638 


IF (IPRINT.NE.4 ) GO TC 6 Cl C 

FMN 

6380 

639 


PRINT 5333* < RT I No* KT ) * KT-1* ISTEP) 

FMN 

6390 

£40 


GC TO 60 1C 

FMN 

640 C 

541 

6203 

CONTINUE 

FMN 

6410 

642 

6303 

CONTINUE 

FMN 

6420 

643 

6403 

CONTINUE 

FMN 

6430 

644 

6503 

CONTINUE 

FMN 

6440 

645 


GO TC 601 0 

FMN 

645C 

646 

C 


FMN 

6460 

647 

c 


FMN 

6470 

848 

c 

KRE ARRAY CONTAINS STANDBY BLOCKS WHOSE RELIABILITY HAS BEEN 

FMN 

646 C 

649 

c 

CALCULATED. 

FMN 

64 30 

650 

c 


FMN 

6500 

651 

6C1C 

DC 6CG4 N-l»i_5 

FMN 

6510 

652 


IF (KRB(N).EQ.O) GO TC 6CG5 

FMN 

6520 

653 

6004 

CONTINUE 

FMN 

6530 

654 

6005 

KRB l M ) -NO 

FMN 

6540 

655 


KRB1-N 

FMN 

6550 

656 

c 


FMN 

656 C 

557 

c 

CHECK IF ALL SENSE BLOCKS ACCOUNTED FOR* IF SO GO TO 8000. 

FMN 

6570 

658 

c 


FMN 

6580 

659 


19=0 

FMN 

6590 

660 


DO 6006 1-1 • HI 

FMN 

660 0 

561 


IF t I0.E3.1 > GO TO 6007 

FMN 

6610 

662 


IF (ISCI ) .EGl.INDEP ) IQ=1 

FMN 

6620 

663 


GC TO 6016 

FMN 

6630 

664 

6007 

ISCI“1)=IS(I l 

FMN 

6640 

665 

6016 

CONTINUE 

FMN 

66 5Q_ 

666 


IF CIStI) .EQ.L) GO TO SDCe 

FMN 

6660 

667 

6006 

CONTINUE 

FMN 

6670 

668 

6008 

NSR-NSR-l 

FMN 

6680 

669 


IF( I.EQ.2) 30 TO 3000 

FMN 

6690 

670 


JJJ- JJJ+1 

FMN 

6700 

671 


IF( JJJ.GE.M4) GO TO 30GC 

FMN 

6710 

672 


GC TC 70 CO 

FMN 

6720 

673 

C 


FMN 

6730 

674 

C 

FINALIZE VARIABLES AND CALCULATE RELIABILITY OF TOTAL SYSTEM. 

FMN 

674C 

675 

c 

SET NOW - 1. 

FMN 

6750 

676 

c 


FMN 

6760 
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677 


acoo 

T-TTOT 


FMN 

6770 

678 



K-NS7EF 


FMN 

678 C 

57 9- 


- - ■ 

N0W=1- -- - ----- - - - - 

- . _ 

F-KN 

67 90 

680 

C 




FMN 

68CC 

631 

C 


DUPLICATE JSAVE INTO ISAVE. 


FMN 

6310 

682 

C 




FMN 

682 C 

583 



CALL QUPTRE C I SAVE t JSltJSAVE* I J5 1 ) 


FMN 

6830 

684 



GO TO 9C0C 


FMN 

6840 

685 

C 




FMN 

6350 

686 

C 


PRINT RESULTS IF RC = 1* SUPPRESS PRINTING OF RO. 

FMN 

68 60 

687 

C 




FMN 

5 3 70 

688 

9CG 1 

PRINT 806 


FMN 

6660 

689 


306 

FORMATC IHi ) 


FMN 

6890 

690 

C 




FMN 

6SCC 

691 



IF t NS8YMX.LE.C) GO TC 700 


FMN 

6910 

692 

C 




FMN 

6920 

693 

C 


PRINT STANDBY BLGCKSC 51- 65) AND THEIR TREES. 


FMN 

6 9 3Q 

6S4 

c 


STORE ALL DIFFERENT STANDEY BLOCKS IN THESE TREES 

IN NRE ARRAY 

. FMN 

6940 

695 

c 




FMN 

6 9 5C 

696 



END FILE ID 


FMN 

6960 

697 



REWIND 1C 


FMN 

6970 

698 



DO 530 N- 1 * N St YM X 


FMN 

6980 

699 



READ(iQ) NS. LJS1 


FMN 

699C 

700 



00 500 I -1 . M2 


FMN 

7C0G 

7C1 



RE A 0(10 CLSAVEC J.I). J-I »M 1 ) 


FMN 

70 IQ 

7C2 


SCO 

continue 


FMN 

702 C 

703 



IF( IPRI.NT.LE. 1) GC TO 52 5 


FMN 

7G3Q 

7C4 



PRINT 510. NS 


FMN 

7C4Q 

7 C 5 


51G 

FORMAT! ' C5TAN0&Y TREE REPLACED 6Y 5L0CK’I3) 


FMN 

7050 

706 



CALL TRFRNT ( S.L S A VE »L JS 1 ♦ C ) 


FMN 

7C&C 

7C7 


525 

CALL TRENUMCLSAVE.LJSl) 


FMN 

7070 

708 



00 520 K-l.NfaNMAX 


FMN 

708 C 

7C9 



NR31-NRDI*! 


FMN 

7G 90 

710 



NRBtNRBl)-NBNUK(K ) 


FMN 

71GC 

711 


520 

CONTINUE 


FMN 

7110 

712 


530 

CONTINUE 


FMN 

712C 

713 

c 




FMN 

7130 

714 

c 


PRINT THE ORIGINAL TREE WITH STANDBY BLGCKS(5i-65)/(Ll-L4). 

FMN 

714C 

715 



IF( IPRXNT.LE. 1) GO TO 700 


FMN 

7150 

716 



PRINT 5 60. LI f L4 


FMN 

716 G 

717 


560 

FORMAT ( * 10R IG IN AL PR08A3ILITY TREE WITH STAND3Y BLOCKS (*I3t 

FMN 

7170 

718 


3 

. ’ THRU* 15. * ).’) 


FMN 

718 C 

719 



CALL TRPRNTtii .ISAVEt US1 . 0 ) 


FMN 

7190 

720 



PRINT 3C6 


FMN 

720C 

721 

c 




FMN 

7210 

722 


700 

PRINT 701 


FMN 

722C 

723 


701 

FORMAK * ACTIVE F/R DORMANT F/ R 

R-INITIAL 

FMN 

7230 

724 


3 

. RELIABILITY’) 


FMN 

724 C 

725 



DO 704 I-1.L4 


FMN 

7250 

726 



IFCI.LE.M1) to TC 700 


FMN 

726 C 

727 



IFC PR031C I ) .LE.O) GO TO 7C4 


FMN 

7270 

726 



PRINT 703. I . PRCBlt I) 


FMN 

7280 

729 


7C3 

FORMAT ( * 3L0CK’I3.44X. 017.71 


FMN 

72 9G 

73C 



GO TO 704 


FMN 

73CC 

7 31 “ 

C~ 

— 





FMN 

.7310.. 

732 

c 


IF IN NR E ARRAY. HAS BEEN REPLACED. THEREFORE NO 1 

RELIABILITY. 

FMN 

7320 

733 

c 




FMN 

7330 

734 


706 

DO 702 J-l . N RD1 


FMN 

734 C 

735 



IFC I. iVE.NRBI J ) ) GO TC 70 2 


FMN 

7350 

736 



IF(RGtl) .GE.i.C » GC TO 711 


FMN 

7 36 C 

737 



PRINT 707. 1$ TUI). TLOCI). ROCI) 


FMN 

7370 

738 



GO TC 704 


FMN 

738 C 

739 


711 

PRINT 709. Ir TLC I) .TLD( I) 


FMN 

7390 

740 



GC TC 704 


FMN 

740 C 

741 


702 

CONTINUE 


FMN 

7410 

742 

c 




FMN 

7420 

743 

c 


DO NOT PRINT SENSE BLOCKS SINCE TLlI)-0» TLDCD-G 

. AND RC< I )-l 

AREFMN 

7430 

744 

c 


NCI RELATED TO THE SENSE SWITCH OPTION PARAMETERS 

AND ARE ONLY 

FMN 

7440 

745 

c 


USED TO SET SENSE BLOCK PROBABILITY = 1.0 IN TREE 

PATHS. 

FMN 

7450 

746 

c 




FMN 

7460 

747 



DC 715 J-l • NN SR 


FMN 

7470 

748 



IFIIISC Jl.EG.I) GC TO 7 C4 


FMN 

7460 

749 


715 

CONTINUE 


FMN 

74 90 

750 

c 




FMN 

75CC 

751 



IFIPROBIC D.LE.Q.) GO TO 704 


FMN 

7510 

752 



lFCRCm-i. ) 7C5 . 7C8 * 7C5 


FMN 

7520 
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753 


7C5 

PRINT 70 7 » I. TLII>* T LD 1 1 > * RO(I). PROBKI) 


FMN 

7530 

754 


7C7 

FORMATS bLO CK * 13 t 2E16.7 *F 12 . 7 t D17 .7 > 



FMN 

7 54 C 

755 



GO TO 704 



FMN 

7550 

756 


706 

PRINT 70S. It TL ( I ) v TLC(I»t PRCBlCIJ 



FMN 

7 56 C 

757 


709 

FORMAT ( * BLOC K # 13 * 2016.7 *12XfD17.7> 



FMN 

7570 

758 


704 

CONTINUE 



FMN 

758 C 

759 

C 





FMN 

7590 

760 

C 


PRINT SENSE SWITCH OPTIONS. 



FMN 

760C 

761 

C 





FMN 

7610 

762 



IF (NNSR.EC.O 1 GO TO 740 



FMN 

7620 

763 



PRINT 720 



FMN 

7630 

764 


720 

F CRH AT ( *0SEN5E SWITCH ACTIVE F/R DORMANT F/R 

PROBABILITY* 

I FMN 

7 64 Q 

765 



DO 730 J-I * NN SR 



FMN 

7650 

766 



NS-I I S ( u ) 



FMN 

766 C 

767 



PRINT 72 5* NS*TLS(NS> *TLDS(NS) *SWPR03tN 

S) 


FMN 

7670 

768 


725 

FORMAT ( * £LG CK *13# 2E16.7 *F 12 . 7 * E 16 . 7 > 



FMN 

768 C 

769 


730 

CONTINUE 



FMN 

7690 

77C 


740 

R-TOT 



FMN 

7700 

771 



K=-L0G1C(1 .-a ) 



FMN 

7710 

772 



K-K+NSIG 



FMN 

7720 

773 



IF ( K.GT.3) K-3 



FMN 

7730 

774 



V ( 12 UVK( K-*2 i 



FMN 

7740 

775 



V ( 14 )-VK ( K ) 



FMN 

7750 

776 



WR IT E ( 6 * V ) TTCTtR 



FMN 

7 76 C 

777 



PRINT 806 



FMN 

7770 

778 

c 





FMN 

7780 

779 

c 


MODIFY PROBABILITY DATA IF LAST = 7* GO 

TO 9BC. 


FMN 

77 90 

780 

c 





FMN 

7800 

781 



IF< LAST. LT. 71 GO TO 998 



FMN 

7810 

782 



IF (LAST. EG. 7 > 60 TO 98C 



FMN 

7820 

783 



IF( LAST, EG. 3) GO TO 999 



FMN 

7830 

784 



IF (LAST .EG. 9 ) GO TO 1C66 



FMN 

7840 

785 



GO TO 933 



FMN 

7850 

786 

c 





FMN 

786 C 

7 87 

c 


INITIALIZE VARIABLES FOR RECALCULATION 1 

WITH NEW PARAMETERS. 

FMN 

7870 

788 

c 





FMN 

788 C 

789 


S8G 

00 990 I-i * Mi 



FMN 

7890 

790 


990 

IS(I>-IISII ) 



FMN 

7900 

791 



NSR-NNSR 



FMN 

7910 

792 



I sue - JSUC 



FMN 

7920 

793 



NOPRNTrQ 



FMN 

7930 

794 



60 TC 997 



FMN 

7940 

735 

c 





FMN 

7950 

796 

c 


ERROR RETURN FROM SUBROUTINES. READ TO 

NEW DIAGRAM 

DATA DECK. 

FMN 

7960 

797 

c 





FMN 

7970 

796 


996 

IF (LAST - 8) 9986*999. 1066 



FMN 

7980 

799 

9998 

READ 270 7 » LAST 



FMN 

7990 

800 

27G7 

FORMAT ( 7 3 X 1 1 1 ) 



FMN 

8C0C 

80 1 



IF ( LAST-8 ) 939 3 •939*1066 



FMN 

8010 

802 

1066 

ENC 



FMN 

8020 
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TREE»FTREE.FLT 


1 


SUBROUTINE FLT(NB.NRSW) 

FLT 

1C 

2 

C 

FAULT TREE AN ALYSIS ~ SU8R OUT INE • " ' " 

FLT- 

- 20- 

3 

C 

UNIVAC 1108 * FORTRAN V VERSION. 

FLT 

3C 

4 

c 

WRITTEN BY PAUL CHELSON. JPL SECTION 153 

FLT 

40 

5 

c 


FLT 

5 C 

6 

c 

THIS SUBROUTINE CONVERTS THE FAULT TREE TO A RELIABILITY 

FLT 

60 

7 

c 

block diagram. 

FLT 

70 

3 

c 


FLT 

80 

9 

c 


FLT 

9C 

10 


C0MM0N/ALLSU3/M1.M2.M3.A 4* M5. M6# M7 • LI *L 3 *L 4 • L5 • L 6 

FLT 

100 

11 


COMM ON /LG SUB/IB 1 15. 2 .50 ) *IIS <5C > • I RB ( 50 1 30 >* HEMP (41 *20 1 t N1 • 

FLT 

lie 

12 


1 IGAVEC 50 #200 ) * ISUC * JS UC . JS 1* JSAVE ( 50*200) • I JS 1 r T OT * PR 08 t 5 5 ) * 

FLT 

120 

13 


2 PRO El 165 J.IPRlNT.NSTDBYdS) .NSBYKX.NOPRNT 

FLT 

13C 

14 


CCMM0N/SU8/IP ATH( 2C.50)* IAND< 49. 12) * IANDZ(25. 1 2 ) 1 1 A N* 8L 3» M 8 

FLT 

140 

15 


CGMM0N/GT/6L1 

FLT 

150 

16 


DIMENSION ISAVE(49) *MZ (49) * A( 8 ) • N3S AVE 1 49) * JPATH (20) 

FLT 

160 

17 


DIMENSION NONACT (50*30) tMATCSW(5Ql 

FLT 

17 C 

18 


EQUIVALENCE ( NONACT (1* 1) *IRB( 1*1) ) * (MATCSU ( 1) * IIS ( 11) 

FLT 

180 

19 


DOUBLE PRECISION FROB* FRO B1 *T0 T 

FLT 

13C 

20 


INTEGER A* BL3 *BLl*GATEf£Q 

FLT 

200 

21 

c 


FLT 

21 0 

22 

c 

IPATH(M*I> HOLOS THE H-TH GATE OF THE I-TH FAULT PATH. (THE I~TH 

FLT 

220 

23 

c 

FAULT PATH IS THE FAULT PATH OF BLOCK I.) 

FLT 

23C 

24 

c 


FLT 

2 4D 

25 


DATA EG/*-*/ 

FLT 

25C 

26 


DATA A/ 0 160 50 5050 5C5 * 0 35 050 50 5D5 05 * 00 60 50 5 0 5C 5 05. 02 40 50 50 50 50 5* 

FLT 

260 

27 


1037C50505C505.03C05C5C5C5C5. C11050505G5C5 .01 40505D505D5/ 

FLT 

27 C 

2a 


2KT RE/ 00505050 50505/ 

FLT 

2 80 

29 


BL3 = KTRE 

FLT 

29 C 

30 


BL1 = KTRE 

FLT 

300 

31 


M 8 - 20 OH A X NC. OF GATES IN A FAULT PATH FOR A BLOCK. 

FLT 

310 

32 


L9 = Ml— 1 dlMAXIMUM NUMBER OF INPUT BLOCKS IN FAULT TREE. 

FLT 

320 

33 


MIC - 12 BM AX NO. OF AN C GATES BEING WORKED ON AT ONE TIME. 

FLT 

33 C 

34 


LL3-L3+1 

FLT 

340 

35 


LL6-M6 ♦ 1 

FLT 

35 0 

36 


L8- M8-1 

FLT 

3 60 

37 

c 


FLT 

370 

33 

c 


FLT 

3 80 

39 

c 


FLT 

39 C 

40 

c 


FLT 

400 

41 

c 


FLT 

410 

42 

c 

INITIALIZE 

FLT 

4 20 

43 

c 


FLT 

430 

44 


898 DO 11 1=1* Ml 

FLT 

4 40 

45 


DO 11 J=1 .2 

FLT 

450 

46 


DO 10 K= 1 • L 3 

FLT 

4 60 

47 


10 IB(K*J*I)=C 

FLT 

470 

48 


DO 9 L-1.M8 

FLT 

4 80 

49 


JP ATH ( L ) - BL 3 

FLT 

490 

5G 


9 IPA TH ( L * I ) -BL 3 

FLT 

500 

51 


11 CONTINUE 

FLT 

510 

52 


DO 12 1= 1 * M7 

FLT 

520 

53 . 


DO 12 J=1.L6 

FLT 

530 

54 



FLT 

- 54C-- 

55 


DO 14 1=1. L9 

FLT 

55 C 

56 


ISA VE 1 I ) -C 

FLT 

5 60 

57 


M Z ( I ) =0 

FLT 

570 

58 


NBSAV£( I ) = 0 

FLT 

5 80 

59 


DO 14 J-1.M1C 

FLT 

590 

60 


14 IAND( I • J ) = 0 

FLT 

6 00 

61 


DO 16 J=1.M1C 

FLT 

610 

62 


I AND ( 1 » J ) - BL3 

FLT 

6 20 

63 


I AND Z ( 1 * J )- BL3 

FLT 

630 

64 


DO 16 K=2 * L 3 

FLT 

6 40 

65 


16 I AND Z ( K • J ) =C 

FLT 

650 

66 


IZ=0 

FLT 

6 60 

67 


MS AVE =0 

FLT 

670 

68 


MBS MAX=0 

FLT 

6 80 

69 


REWIND 10 

FLT 

690 

70 

c 

INITIALISE VARIABLES USED IN MODIFYING TREES WITH DORMANCY. 

FLT 

7 00 

71 


DO 517 J = 1.50 

FLT 

71C 

72 


MATCSW(J) = 0 

FLT 

7 2C 

73 


DO 517 K = 1 .30 

FLT 

73 C 

74 


NONACT ( J *K ) = 0 

FLT 

7 40 

75 


517 CONTINUE 

FLT 

750 
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76 

C 



FLT 

7 6C 

77 

C 



FLT 

77 C 

78 

c 


PHASE 2 

FLT 

7 80 

79 

c 


READ FAULT TREE* WRITE CN SCRATCH TAPE. STORE BLOCK NUMBERS IN 

FL T 

79C 

3D 

c 


NBSAVE ARRAY. N3KMAX - MAX NUMBER BLOCKS IN DIAGRAM. 

FLT 

8 00 

81 

c 


NB-BGTTCM BL OCK »K -PAT M KO..U=A.I.OR OtI=A»ItOR C IDENTIFICATION NOFLT 

81 C 

82 

c 


RANK N3SAVE. LOWEST TO HIGHEST. ELIMINATING EQUAL BLOCK NUMBERS. 

FLT 

8 20 

83 

c 



FLT 

830 

84 

c 



FLT 

8 40 

85 


2C 

READC5*1CC*ERR=1C5*END=1C5) N 6 t JF AT Ht LAS T 

FLT 

85C 

86 


100 

F0RMATCl2.20A3.17X.il) 

FLT 

8 60 

87 



NBSMAX=NBSMAX+1 

FLT 

870 

88 



NBSAVEl N8SMAX >-N3 

FLT 

8 80 

89 



GO TO 514 

FLT 

89 0 

9G 


105 

READ(C.i06»ERR=104.END=104) I.J 

FLT 

900 

91 


106 

F0RMATC7X .11 .6X.I1 ) 

FLT 

91 C 

92 



IFC CI.GT.4 J.OR.IJ.GT.a )> GO TO 1C4 

FLT 

920 

93 



IFCJ.EG.O) J=3 

FLT 

93C 

94 



NSI G- J 

FLT 

940 

95 



IFRINTzI 

FLT 

950 

96 



PRINT 109.IPRINT.NSIG 

FLT 

960 

37 


109 

FORM ATI *DI PRINT OPTION SET AT»I4.* NSIG OPTION SET AT*I4/1H1) 

FLT 

37 C 

98 



GO TO 20 

FLT 

980 

99 


104 

PRINT 108 

FLT 

990 

100 


108 

FORMATI *QERRGR** ATTEMPT TO SPECIFY IPRINT OR NSIG. BUT SPECIFIED 

FLT 

1CGC 

1D1 


10 UT OF RANGE.*/* CARD IGNORED*) 

FLT 

1C1C 

1C 2 



GO TO 20 

FLT 

1020 

1C3 


514 

WRITE CIO .NB.JPATH.LAST 

FLT 

1030 

104 



IFCLAST.LT.l) GO TO 20 

FLT 

1040 

105 



END FILE 1C 

FLT 

1C5G 

108 



REWIND 10 

FLT 

1C60 

1C7 

c 



FLT 

1C7G 

1C8 

c 



FLT 

1C 80 

1C9 

c 


NBSAVE ORIGINALLY HOLDS BLOCK NCS USED BY THE FAULT TREE. IT IS 

TFLT 

1090 

110 

c 


NBSAVE ORIGINALLY HOLDS BLOCK NOS USED 3Y THE FAULT TREE. IT IS 

FLT 

1100 

111 

c 


THEN FILLED CUT WITH NOS. NOT USED. THUS GIVING A LIST OF NUMBERS 

FLT 

1110 

112 

c 


THAT CAN BE ASSIGNED TO EQUIVALENT BLOCKS. 

FLT 

1120 

113 

c 



FLT 

113 C 

114 



IFC IPRINT. NE. 4) GO TO 205 

FLT 

1140 

115 

c 


PRINT NBSAVE TO CHECK 

FLT 

115 G 

116 

c 



FLT 

1160 

117 



PRINT 201.NBSAVE 

FLT 

117 C 

118 


201 

FORMAT! *0NBSAVE*25I3/1H » SX. 2513) 

FLT 

11 8C 

119 

c 



FLT 

119 C 

120 

c 



FLT 

1200 

121 

c 


NBKMAX = MAX NUMBER OF ELCCKS IN DIAGRAM. 

FLT 

121 C 

122 

c 



FLT 

1220 

123 


205 

NBKMAX=NBSMAX 

FLT 

1230 

124 

c 



FLT 

1240 

125 

c 


RANK NBSAVE 1 BLOCK NUMBERS) ELIMINATING EQUALS. 

FLT 

125C 

126 

c 


JCOUNT = NUMBER OF EQUIVALENT BLOCKS. 

FLT 

1260 

127 

c 



FLT 

127 C 

123 



JCOUNT-O 

FLT 

1280 

129 



IEND-NBKMAX-1 

FLT 

123 C 

130 



DO 250 I-l.IEND 

FLT 

1300 

131 



JBEG -I+i 

FLT 

1310 

132 



DO 240 J=UBEG»N8KMAX 

FLT 

1320 

133 



IF CNBSAVEC J).EQ.C) GO TC 25C 

FLT 

133 C 

134 



IFC NBSAVEC I > . EQ.N3S AV£ < J )) GO TO 200 

FLT 

1340 

135 



IF CNBSAVECI ) .LT.NBSAVEC J) ) GC TO 24C 

FLT 

1350 

136 



JH0LD-N33AVEC I> 

FLT 

1360 

137 



NBSAVECI)-NBSAVE(J) 

FLT 

137 C 

138 



NBSAVEC J)-J HOLD 

FLT 

138C 

139 

140 

c 


GO TC 240 

FLT 
- F LT“ 

1390 

T4C0 

141 

c 


“CHECK IF DUPLICATE ALREADY IN 1TEMP . IF NOT STORE. 

FLT 

1410 

142 

c 



FLT 

1420 

143 


200 

DO 2 1C L-1.M7 

FLT 

1430 

144 



IFC ITEMPC2.L) .£Q.0> GO TO 215 

FLT 

1440 

145 



IF CITEHPC2»L).EQ.NBSAVECI1> GO TO 225 

FLT 

1450 

146 


210 

CONTINUE 

FLT 

1460 

147 



PRINT 211 

FLT 

1470 

148 


211 

FORMAT! *1ITEMP IS FULL-S TAT ME NT NC. 211 OF FLT SUBROUTINE.*/* THERFLT 

14 80 

149 



IE ARE MORE THAN 20 SETS OF EQUIVALENT BLOCKS. •) 

FLT 

1490 

150 



GO TO 9999 

FLT 

1500 

151 


215 

ITEHPC2.L1-NBSAVEC I> 

FLT 

1 51G 
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152 

C 







FLT 

1520 

153 

_ c 


ELIMINATE DUPLICATE ELO CK NUMBERS AND KEEP COUNT 

• 

FLT 

1 53 C 

154 

c 


SHIFT N3SAVE ARRAY TO LEFT. ’ 

— - 

- 

- - 

- 

FLT 

154C 

155 

c 







FLT 

155 C 

156 


225 

00 230 K-J.IEND 





FLT 

1560 

157 



N6SAVE (K ) -NBSAVECK41 ) 





FLT 

157 C 

158 


2 30 

CONTINUE 





FLT 

1580 

153 



N6SAVEINBKMAXJ-C 





FLT 

159 G 

1 50 



JCOUNT- JCOUNT *1 





FLT 

1600 

161 


240 

CONTINUE 





FLT 

1 61 C 

162 


250 

CONTINUE 





FLT 

1620 

163 

c 







FLT 

1 63 C 

164 

c 


N8EG - END OF ORIGINAL BLOCK NUMBERS IN 

NBSAVE 

ARRAY. 

FLT 

1640 

165 

c 


NEND = END OF ORIGINAL ELCCK NUMBERS IN 

NESA VE 

ARRAY . 

FLT 

1 65 C 

166 

c 







FLT 

1660 

167 



NBEG-NBKM AX- JCOUNT 





FLT 

1670 

168 



NEND-NBEG 





FLT 

168C 

163 

c 







FLT 

169 C 

170 

c 


PUT NEW 8L0CK NUMBERS INTO LAST PART OF 

NBSAVE 

ARRAY. 

FLT 

1700 

171 

c 







FLT 

1 71 C 

172 

c 







FLT 

172C 

173 

c 


PRINT NBSAVE TO SEE IF RANKE C 





FLT 

173 C 

174 



IF! IPRINT.NE. 4) GO TO 271 





FLT 

1740 

175 

c 







FLT 

1 75 C 

176 



PRINT 201 • N8S A Vl 





FLT 

1760 

177 

c 







FLT 

1 7 7 C 

178 


271 

K-l 





FLT 

1780 

179 



N-NBEG 





FLT 

1 79 C 

180 



DO 280 J-l » L9 





FLT 

1 8CC 

181 



IF !NBSAVE!K)*EG.Jl GO TC 270 





FLT 

1810 

182 



N-N«-l 





FLT 

1820 

183 



IFIN.GT.L91 GO TO 26C 





FLT 

183 C 

184 



NBSAVE ( N )zj 





FLT 

1840 

185 



GO TO 280 





FLT 

185 C 

186 


270 

K- K ♦ 1 





FLT 

I860 

187 



IF (K.GT.NEND1 GO TO 265 





FLT 

187 C 

183 


280 

CONTINUE 





FLT 

1 8 8C 

183 



GO TC 29C 





FLT 

189 C 

190 


260 

PRINT 261# N 





FLT 

1900 

1S1 


261 

FORMAT! *1N IS GREATER T hAN 5C. 

N - * 1 4 / • 

ERROR 

AT 

STATEMENT 

NO. 2ECFLT 

191 C 

192 



1 OF FLT SU3R0UTINE.* i 





FLT 

1920 

193 



GO TO 9999 





FLT 

193 C 

194 


265 

LBE G- J + l 





FLT 

1940 

195 



DO 275 L-LBE Gt L 9 





FLT 

1950 

196 



R=N*1 





FLT 

1960 

197 



IFCN.GT.L9) GO TO 260 





FLT 

197 C 

198 



NBSAVE ( N 1-L 





FLT 

1930 

199 


275 

CONTINUE 





FLT 

199 C 

200 

c 







FLT 

2C0G 

201 

c 


PRINT NBSAVE TO CHECK R /NK AND 

NEW NE ASSIGNED 

• 


FL T 

201 C 

202 

c 







FLT 

2C2G 

203 



IF (IPRINT.NE.4) GO TO 290 





FLT 

2 03 C 

204 



PRINT 201 v NBSAVE 





FLT 

2040 

205 — 

c- 

_ 







FLT 

205 C 

20 6 

c 


REREAO DATA AND CHECK IF NB EQUAL I TEMP. 


" - 



FLT- 

206G- 

2C7 

c 







FLT 

2070 

208 


290 

READ! 10 ) NK t JPATH# L AS T 





FLT 

2080 

209 



DO 300 L=1#M7 





FLT 

2090 

210 



IF! ITEMPt 2tL) .EQ.OJ GO TO 30C 





FLT 

2100 

211 



IF (I TEMP!2»L ) • E Q • NK ) GO TC 31C 





FLT 

211 C 

212 


300 

CONTINUE 





FLT 

2120 

213 



GO TC 350 





FLT 

213 G 

214 

c 







FLT 

2140 

215 

c 


CHECK IF FIRST TIME R E A C IN . 

IF NOT AS 

SIGN NEW 

N6. 

FLT 

2150 

216 

c 







FLT 

2162 

217 


210 

ITEMPll *L l-ITE.MPil*L Ul 





FLT 

2170 

218 



IF! ITEMPi 1#L) .LE.ll GO TO 3 50 





FLT 

2180 

219 



NEEG -NBEG+l 





FLT 

2 19 C 

220 



NS-NBSAVE! N8EG) 





FLT 

2200 

221 



N- ITEMP ! 1 »L 1 4l 





FLT 

2210 

222 



ITEMP! N»L)-NB 





FLT 

222C 

223 



GO TO 355 





FLT 

223 G 

224 

c 







FLT 

2240 

225 

c 


STORE IPATH WITH NEW NB IF NEEDED. 




FLT 

2250 

226 

c 







FLT 

226C 

227 


350 

NB=NK 





FLT 

22 7 G 
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228 


355 

DO 360 1=1 * M8 


FLT 

2280 

229 


360 

IP ATHtI*NB) = JPATH! 11 


FLT 

2 29C 

230 



IF ( LAST # LT * 1 ) 60 TO 290 


FLT 

2 3CC 

231 

C 




FLT 

2 31C 

232 

C 


REINITIALIZE ITEMP 


FLT 

2 32C 

233 

C 




FLT 

2 33 C 

234 



DO 370 L-1.M7 


FLT 

23 4G 

235 


370 

ITEMF<1*L)=0 


FLT 

235 C 

236 

c 




FLT 

2360 

237 

c 


PRINT EQUIVALENT BLOCKS. N1 - NUMBER OF CF EQUIVALENT BLOCK 

SETS 

*FL T 

2 27 C 

233 

c 




FLT 

2380 

239 



PRINT 371 


FLT 

239 C 

240 


371 

FORMAT! *OEQUI VALENT BLOCKS*) 


FLT 

2400 

241 



DO 380 I-1.M7 


FLT 

2 41 C 

242 



IF! ITEMPI2.I) .EQ.O) GO TO 390 


FLT 

242C 

243 



DO 373 J=2tLL6 


FLT 

243 C 

244 



IF! ITEMPC Jtll *E3.0) GO TO 374 


FLT 

2440 

245 


373 

CONTINUE 


FLT 

245 C 

246 


374 

K= J— 1 


FLT 

2460 

247 



PRINT 372 »I T EHP ! 2 • I ) * ! ! £Q *ITEHP!J*I) 1 * J = 3 *K 1 


FLT 

247 C 

248 


372 

FORMAT ( * 0* 20 ( 13* Al) ) 


FLT 

2480 

249 



DC 380 K=2.LL6 


FLT 

249 C 

250 



ITEMP! K*M6*I)=-ITEMP!K*I ) 


FLT 

2500 

251 


380 

CONTINUE 


FLT 

251 C 

252 



Nl = M7 


FLT 

2520 

253 



GO TO 21 


FLT 

253 C 

254 


390 

Nl- I— 1 


FLT 

2540 

255 



IF ! Nl* E Q .0 1 PRINT 391 


FLT 

255 C 

256 


391 

FORMAT ( *0 NONE USED* ) 


FLT 

2560 

257 

c 




FLT 

257 C 

258 

c 


PRINT FAULT TREE PATHS AS A CHECK 


FLT 

2580 

259 

c 




FLT 

259 C 

260 


21 

PRINT 633 


FLT 

2600 

261 


633 

F0RMAT11H116HFAULT TREE PATHS) 


FL T 

281 C 

262 



DO 90 I- 1 * Ml 


FLT 

2620 

263 



IF CIPATH! 1*1 1-BL3 191 *9C .91 


FLT 

2630 

264 


91 

PRINT 634. It ( IPATH(K.I )♦ K-1.M8) 


FLT 

2640 

265 


634 

FORM AT! 1H05HBL0CKI3.3X* 20 IA3 •1X1) 


FLT 

265 C 

266 


90 

CONTINUE 


FLT 

2660 

267 

c 




FLT 

267 C 

268 

c 


PRINT ITEMP AS A CHECK IF IPR INT IS 4 


FLT 

26 80 

269 

c 




FLT 

269 0 

270 



IF! IPRINT.NE* 4) GO TO 410 


FLT 

2700 

271 



DC 400 I=1*M7 


FLT 

271 C 

272 



PRINT 393. I* IITEMPCJ*!). J=1.L6) 


FLT 

2720 

273 


393 

FORMAT! *CI TEMP ( * 13 * * ) = *2113/* *20131 


FLT 

273 C 

274 


4 00 

CONTINUE 


FLT 

2740 

275 

c 




FLT 

275 C 

276 

c 




FLT 

2760 

277 

c 


PHASE 3 


FLT 

277C 

278 

c 




FLT 

2780 

279 

c 


THE SUBROUTINE MODIFY MODIFIES THE FAULT PATHS WHERE NECESSARY 

FLT 

279C 

280 

c 




FLT 

2800 

281 


410 

CALL MODIFY! NO NACT tMATCSW »A*IPRINT*NRSW) 


FLT 

2810 

282 

c 




FLT 

2820 

283 

c 


THE MODIFIED FAULT TREE HAS ELIMINATED ALL SERIES *ANC* AND 

♦ OR* 

FLT 

2 83 C 

284 

c 


GATES* AND HAS REOUCED ‘STANDBY* GATES TO * AND * AND *0R* GATES 

FLT 

2840 

285 

c 


! BUT WHICH ARE DENOTED 0 AND S RESPECTIVELY! 


FLT 

285 C 

288 

c 




FLT 

2860 

287 


578 

IFlIPRINT.EG.O) GO TO 5 C8 


FLT 

287 C 

288 



PRINT 577 


FLT 

2880 

289 


577 

FORMAT (1H1*5X* * MODIFIES FAULT TREE PATHS*) 


FLT 

289 C 

290 



DO 505 1=1 9 Ml 


FLT 

2900 

291 



IF ! IPATH(1*I) -BL3) 50E*5C5*5C6 



-FLT 

-291-G 

2.92 

5Q6 — 

-PRINT— 5Q7*“I-tl-rP ATH TKTI )VK= 1 * M 8 ) 


FLT 

2920 

293 


507 

FORM AT! 1HC5HEL0CKI3 *3X* 2C !A3 *1X ) ) 


FLT 

293 C 

294 


505 

CONTINUE 


FLT 

294C 

295 

C 




FLT 

295 C 

296 

C 


PHASE 4 


FLT 

2960 

297 

C 


CONVERT CONDENSED PATHS TO RELIABILITY BLOCK DIAGRAM 


FLT 

2 970 

298 

c 




FLT 

2980 

299 


508 

I K =0 


FLT 

299C 

300 



I AN=0 


FLT 

3000 

301 



KS = D 


FLT 

3010 

302 

c 




FLT 

3C20 

303 

c 


IK ITERATES FAULT PATHS 


FLT 

3030 
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c 


304 

305 

306 

307 
303 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

332 

333 

334 

335 

336 

337 

338 

339 

340 

341 


C 

C 

C 

C 


C 

C 

C 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


40 IK-IK + l ------- ------ _ _ 

HK-0 

IF IIK-M1 ) 1 .82 .82 

MK ITERATES THE GATE IN THE FAULT PATH 
1 HK-MK+1 

IF { MK-M8 )170t 170.40 
170 IF CIPATHIMK. IK > -BL3 > 4 1 . 40 .41 

IS IT X 

X MEANS THAT THIS GATE HAS ALREADY BEEN TREATED 

41 IF { IPATH( MK * IK ) — A ( 2 ) >56.1.56 

56 J = IP ATH ( MK .IK ) 

57 MSAVE^l 

42.43 LOOP THROUGH ALL OTHER FAULT PATHS TO SEE IF GATE J IS IN 
ANY OTHER PATH. IF YES. LOOP CUT TO 61. 

DO 42 1=1. L9 

IS PATH TO BE CHECKED BLANK * IF YES. IGNORE 
IF C IP ATH (1. 1 J-BL3J58.42 .58 

IS THIS THE PATH CURRENTLY BEING WORKED ON (IK) * IF YES t IGNORE 

58 IF(I-IK>59. 42.59 

59 DC 43 N=1 . MS 

IF (IFATHCN.I J-J )43.6l*42 
43 CONTINUE 
GO TO 42 

IS IT CLEARED 

IS GATE FOUND THE FIRST GATE IN THE FAULT PATH. IF NOT. HAS THE 
PREVIOUS GATE IN THE FALLT PATH BEEN TREATED * I.E. IS THE 


FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 

FLT 


342 

C 

PREVIOUS GATE = X ♦ 




FLT 

343 

C 





FLT 

344 

61 

IF ( N- 1)62*63*62 




FLT 

345 

62 

IF ( IP ATH ( N-l *II-A( 2 1 )64 .6 3.6 4 




FLT 

346 

C 





FLT 

347 

C 

NOT CLEARED IS 64 ROUTE 




FLT 

348 

C 





FLT 

349 

64 

IF(IZ)172. 518.172 




FLT 

350 

17 2 

IF ( G ATE ( J ) — A I 7 ) ) 174.176.174 




FLT 

351 

174 

IF (GATEC J)-A(4 > ) 173.175.173 




FLT 

352 

173 

IF (GATE(J) - A ( 6 ) ) 176.175.176 




FLT 

353 

176 

1Z=0 




FLT 

354 

518 

MK=0 




FLT 

355 


1K = I 




FLT 

356 


GO TO 1 




FLT 

357 

~c 




— 


_FL T 

358 

c 

CHECK TO SEE IF OR GATE HAS BEEN 

1 HIT 

BEFORE 

FLT 

359 

C 





FLT 

360 

175 

DO 140 K=1 # Ml G 




FLT 

361 


IF IIANDZC I.K )-J)14C»I78 .140 




FLT 

362 

140 

CONTINUE 




FLT 

363 


DO 146 K = 1 1 MIC 




FLT 

364 


IF ( IANDZI 1.KJ-BL3) 146. 177.146 




FLT 

365 

146 

CONTINUE 




FLT 

366 

C 





FLT 

367 

C 

IANZ IS FULL. IMPLYING THAT MORE 

THAN 

12 

OR GATES ARE INTER- 

FLT 

368 

C 

RELATEO. J IS THE OR GATE 8EING 

WORKED 

ON. 

FLT 

369 

C 





FLT 

370 

5 

PRINT 15. J 




FLT 

371 

15 

F ORM AT ( • 1 ERROR ♦♦ THERE ARE MORE 

THAN 

12 

OR GATES INTERRELATED.*/ 

FLT 

372 


L* THE GATE CURRENTLY BEING WORKED 

ON 

IS 

• A 3/ * ERROR IS AT STATEMENFLT 

373 

2T NO. 5 OF FLT SUBROUTINE ) 




FLT 

374 


GO TO 9999 




FLT 

375 

C 





FLT 

376 

177 

IAN— K 




FLT 

377 


I AND Z { 1 .1 AN ) =J 




FLT 

378 

C 





FLT 

379 

C 

PUT BLOCKS WITH TEMPORARY OUTPUT 

CF Z 

INTO I AN02 ( ..IAN) 

FLT 


3040 
305 C 
30 60 
307 C 
30 80 
3C9C 
31G0 
3110 
3120 
313 G 
3140 
315 C 
3160 
3170 
3180 
319 D 
3200 
3210 
32 20 
323C 
3240 
3250 
3260 
327C 
32 80 
3290 
33C0 
331C 
3320 
333 C 
3340 
335 C 
3360 
3370 
3380 
333 C 
3 4 CO 
3410 
34 20 
3430 
34 40 
345 C 
3460 
347C 
3480 
349 C 
3500 
3510 
3520 
353 C 
3540 
355 C 
3560 
357 C 
3580 
359 C 
3600 
361 C 
3620 
363 C 
3640 
365 C 
3660 
3670 
3680 
369 0 
3700 
371 C 
37 20 
373C 
3740 
375 C 
3760 
377 C 
37 80 
379C 
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380 

C 







FLT 

3QCG 

381 


4 

IFIIZ.LT.L3J 60 TO 1004 





FL1 

381 C 

382 



STOP 1004 





FLT 

3820 

383 

1004 

DO 141 K=I*IZ 





FLT 

383 C 

384 


141 

IANDZ<K+i.IAN)=MZ<K) 





FLT 

33 40 

385 



60 TO 176 





FLT 

385 C 

386 


178 

IAN-K 





FLT 

38 6C 

387 


7 

CALL SU E7 ( I B ) 





FLT 

387 C 

388 


149 

I AND! 1 . I AN ) -BL 3 





FLT 

3880 

389 



DO 144 K = 2 »L3 





FLT 

389 C 

390 



IANDZ(K.IAN)=0 





FLT 

3 90Q 

391 


144 

I AND l K »I AN ) = Q 





FLT 

391C 

392 



DO 145 K - L L 3 » L9 





FLT 

3920 

393 


145 

I AND (KtlAN)-C 





FLT 

393 C 

394 



GO TO 4 





FLT 

3940 

395 

C 







FLT 

395 C 

396 

c 


CLEARED IS 63 ROUTE 





FLT 

3960 

397 

c 







FLT 

397 C 

398 


63 

MSAVE=MSAVE*1 





FLT 

3980 

399 

c 







FLT 

399C 

400 

c 


ISAVE SAVES PATH NUMBERS OF THE 

PATHS 

THAT 

CONTAIN GATE 

: j 

FLT 

4000 

401 

c 


MSAVE IS THE NUMBER OF ELOCKS STORED 

IN ISAVE 


FLT 

4D1C 

402 

c 







FLT 

4020 

403 



ISAVE(MSAVE)-! 





FLT 

403 0 

404 


42 

CONTINUE 





FLT 

4040 

405 

c 







FLT 

4C5C 

4C6 

c 


RECORD ORIGINAL PATH 





FLT 

40 60 

407 

c 







FLT 

4C7C 

408 



ISAVEC1)=ZK 





FLT 

4C8C 

409 

c 







FLT 

4C9C 

41C 

c 


SET THE FIRST NON-X-ED GATE IN 

EACH OF THE 

PATHS SOREO 

IN ISAVE 

FLT 

4100 

411 

c 


TO X. THIS INDICATES THAT 

THESE 

GATES 

ARE WORKED 

ON. 

FLT 

4110 

412 

c 







FLT 

4120 

413 



DC 44 1 = 1. MSAVE 





FLT 

413 C 

414 



L=ISAVE( I) 





FLT 

4140 

415 



DC 45 N=1.M8 





FLT 

4150 

416 



IFC IP ATH < N»LJ-AC2>>55. 45. 65 





FLT 

4160 

417 


45 

CONTINUE 





FLT 

417C 

413 


65 

IPA TH ( Nt L )= A < 2 ) 

* 




FLT 

4180 

419 


44 

CONTINUE 





FLT 

419 C 

420 

c 







FLT 

4200 

421 

c 


RECORD ISAVE IN I AND FO F ANY OR 

GATES BEING WORKED ON 


FLT 

421 C 

422 

c 







FLT 

422C 

423 



DO 15C 1=1. MIC 





FLT 

423C 

424 



IF ( I ANDZ( 1 »I J—BL 3)182*1 50 *13 2 





FLT 

4240 

425 


182 

IFlIANDZd.I J-JJ79C.150 .790 





FLT 

425 0 

426 


790 

DO 151 K-2 »L9 





FLT 

4260 

427 



IF ( I ANDlKtini51 .183.151 





FLT 

4270 

428 


151 

CONTINUE 





FLT 

4280 

429 


183 

IAND(l*I)=IANDZtltI> 





FLT 

429 C 

430 



IL1 - MSAVE ♦ K - 1 





FLT 

4300 

431 



DC 152 N=K.lLl 





FLT 

4310 

432 



L2-N-K+1 





FLT 

4320 

433 


152 

IANDCN.I»=ISAVE<L2> 





FLT 

4330 

434 


150 

CONTINUE 





FLT 

4340 

435 

c 







FLT 

4350 

436 

c 


NOTE IN THE FOLLOWING*. MZ HOLDS BLOCKS OF 

THE DIAGRAM 

THAT HAVE 

FLT 

4360 

437 

c 


A TEMPORARY OUTPUT OF Z« I .E . OUTPUT 

BLOCK 

NOT YET DETERMINED. 

FLT 

437 C 

438 

c 


IZ IS QTY IN MZ. 





FLT 

4380 

439 

c 

TOST FOR A • 0 OR S.O BY WHETHER LT 

OR GT 

*1* . 



FLT 

439C 

440 

c 







FLT 

4400 

441 



IF ( GATE! J) - Atl>» 67.67.66 





FLT 

4410 

442 

c 







FLT 

4420 


443 C AND GA TES USE 67 RO UTE* FLT 4 43C 


444 

c — 

— 



FLT 

4440 

445 

67 


DC 50 1 = 1. MSAVE 


FLT 

4450 

446 



L=ISA VE( I > 


FLT 

4460 

447 



IB (2 .2 «L ) =A ( 5) 


FLT 

4470 

446 


50 

MZ< I ) =L 


FLT 

4480 

449 



I Z=HS A VE 


FLT 

4490 

450 



GO TO 1 


FLT 

4500 

451 

C 




FLT 

4510 

452 

C 


HAS THIS OR GATE APPEARED IN IANDZ 


FLT 

4520 

453 

c 




FLT 

453C 

454 


66 

DO 153 1=1. M10 


FLT 

4540 

455 



IF! IANDZ lit I l-J 1153 .184 .153 


FLT 

4550 


18 


JP Li Technical Memorandum 33-512 



456 


153 

CONTINUE 

FLT 

456C 



- - - , - -457 


185 

CONTINUE. 

flt 

4 57 C 



453 



IFCIZ )7Q *1 88# 7C 

FLT 

4 580 



459 

c 



fl t 

4 59 C 



460 

c 


PUT BLOCKS IN SERIES. PUT Z IN OUTPUT OF LAST. 

FLT 

4600 



461 

c 



FLT 

4 61 C 



462 


188 

L^ISAVEi 1 J 

FLT 

462C 



463 



1612 t2*L )-ISAVE(2) 

FLT 

4 63 C 



464 



IF(MSAV£>2)73t73f72 

FLT 

4640 



465 


72 

N- MS AVE “ 1 

FLT 

465 C 



466 




FLT 

466G 



467 



DO 47 1-2 *N 

FLT 

4 67 C 



468 



L^ISAVEt 1-1) 

FLT 

468C 



469 



M- IS AVE 1 1 +1 ) 

FLT 

4 69 C 



470 



LM- ISAV£( I ) 

FLT 

4 7CC 



471 



1612 *1*LM)-L 

FLT 

4 710 



472 


47 

18 ( 2 * 2 # LM 1 - M 

FLT 

4720 



473 


73 

LM-ISAVE(MSAVE) 

FLT 

4 73 C 



474 



L-ISAV£( MSAVE-1) 

flt 

4740 



475 



IE 12 *1 * L M >-L 

FL 1 

4 75 C 



476 



I3( 2t 2*LM)-A< 5) 

FLT 

4760 



477 



IZ=1 

FLT 

477 C 



478 



MZ< D-LM 

FLT 

4780 



479 



MK-1 

flt 

4 79 C 



480 



IK-LM 

flt 

4800 



481 

c 



flt 

4810 



482 

c 


blank out paths with outputs 

FLT 

4 82C 



483 

c 



FLT 

4 83 C 



484 


793 

IL1-MSA VE— 1 

FLT 

4840 



485 



DC 46 1-1 *1 L 1 

FLT 

485 C 



486 



LcISAVEi I) 

FLT 

4860 



487 



DO 46 H=1*M8 

flt 

4 87 C 



488 


46 

IPATH( M t L 1- 3L 3 

FLT 

4 8 8C 



489 



60 TO 1 

flt 

4890 



490 


184 

IAN-I 

FLT 

4 9 CC 



491 



IF (I AND (1*1 1 -BL 3 1186. 16 7 *166 

FLT 

491 C 



492 


186 

CALL SU37< IB) 

FLT 

4920 



493 

c 



FLT 

493 C 



494 

c 


REMOVE IANOZ FROM ISAVE 

FLT 

4940 



495 

c 



FLT 

4 95 C 



496 


187 

00 154 Ir2*L 3 

FLT 

4960 



497 



IF ( IAN DZ( I* I AN 1)51*15 4* 51 

FLT 

4 97 C 



498 


51 

DO 155 N- 1 * MS A V£ 

FLT 

493C 



499 



IF(IANCZ( It I AN 1 -ISAVE (N 1)155 *1SL* 155 

FLT 

4990 



SCO 


155 

CONTINUE 

FLT 

5CC0 



501 



60 TG 154 

FLT 

5C1C 



502 


190 

IF ( N-MSAVE 1 18 9 f 19 1# 189 

flt 

502C 



503 


191 

ISAVE(MSAVE)-C 

FLT 

5C3C 



504 



MSA VE-MSAVE-1 

FLT 

5C4C 



505 



GO TC 154 

flt 

5C5C 



506 


189 

IL1 -MSAVE-1 

FLT 

5060 



507 



DO 156 K- N f I Li 

FLT 

5070 



508 


156 

ISAVE( M)-ISAVE(M+1) 

FLT 

50 8G 



509 



-MS AVE-MS AVE“ 1 

FLT 

5C9C 



510 


154 

CONTINUE ' ' 

FLT 

5100 

— 

— 

511 

c 

BLANK 1 AN 0 Z * 1 AND ( » 1 1 AN 1 FERE. 

FLT 

511C 



512 



iAND(iflAN) = BL3 

FLT 

5120 



513 



I AND 2 ( 1 *1 AN ) = BL3 

FLT 

513C 



514 



00 867 K - 2#L3 

FLT 

5140 



515 



I AND Z ( K • 1 AN 1 =■ 0 

FLT 

515 C 



516 



lAND(KflAN) - C 

FLT 

5160 



517 


867 

CONTINUE 

FLT 

5170 



518 



DO 868 K = LL 3 *L9 

FLT 

5180 



519 



IAND CKflAN) - 0 

FLT 

5190 



520 


868 

CONTINUE 

FLT 

520C 



521 


70 

N = Q 

FLT 

521 C 



522 



IZ=0 

flt 

522C 



523 

c 



FLT 

5230 



524 

c 


Eliminate blocks with output from isave 

FLT 

5240 



525 

c 



FLT 

525 C 



526 


74 

NzN*1 

FLT 

52 60 



527 



IF (MSAVE-N)75*76t76 

FLT 

527 C 



528 


76 

L- ISAVE ( N ) 

FLT 

52 80 



529 



IF I IB (2 *2 fL ) 177*74*77 

FLT 

52SC 



530 


77 

IZ=IZ*1 

FLT 

53CC 



531 



MZ(IZI-L 

FLT 

531 C 
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532 



L-N-l 

533 


78 

l~l*i 

534 



ISAVE<L) = ISAV£(L*H 

535 



IFCL«1-MSAVE>78 .79*79 

536 


79 

MSAVE-MSAVE-1 

537 



N = N~1 

538 



IF < MSAVE)74#1 .74 

539 


75 

L-ISAVEC1) 

540 

C 



541 

C 


HAKE CHAIN 

542 

C 



543 



DO 48 1=1. IZ 

544 



M= M 2 ( I > 

545 

C 



546 

c 


BLANK OUT PATHS WITH NUMERICA 

547 

c 



548 



DO SCO K=1*M8 

549 


800 

IPATH( KtMU6l3 

550 



N = I*1 

551 



18 f 2 » 2 • M ) = L 

552 


48 

IB<N.i.L»=M 

553 



IFC MSAVE-1 )81 .81.80 

554 


80 

DO 49 1 = 2 *MS AVE 

555 



L=ISAVE(I-1> 

556 



M=ISAVEtr l 

557 



ISI 2.2»L)=M 

558 


49 

1612 . 1 »M )=L 

559 


81 

L= I SAVE ( MS AVE 1 

560 



IB(2»2.L)=A(5k 

561 



12-1 

562 



MZ (1 1 = L 

563 



IK=L 

564 



MK = 0 

565 



IFCMSAVE-1U. 1.793 

566 

c 



567 

c 


SET UP Ml AS TOP BLOCK 

568 

c 



569 


82 

DO 52 1=1.12 

570 



L=M2CI> 

571 



I1 = 3>1 

572 



1612 .2 »L ) =M1 

573 



NS = Ml 

574 


52 

IB(I1.1.M1>=L 

575 

c 



576 

5999 

continue 

577 



REWIND 10 

578 

10CC0 

RETURN 

579 



END 


FLT 532C 
FLT 533 C 
FLT 534C 
FLT 535 C 
FLT 5360 
FLT 537 C 
FLT 538C 
FLT 539 C 
FLT 5400 
FLT 541 C 
FLT 54 20 
FLT 543C 
FLT 5440 
FLT 54 5 C 

OUTPUT FLT 546C 

FLT 547C 
FLT 548 C 
FLT 54 90 
FLT 550 C 
FLT 5510 
FLT 552 C 
FLT 5530 
FLT 554 C 
FLT 5550 
FLT 558 C 
FLT 5570 
FLT 558 C 
FLT 5590 
FLT 56CC 
FLT 5610 
FLT 562 C 
FLT 5630 
FLT 564 C 
FLT 5650 
FLT 566 C 
FLT 5670 
FLT 568 C 
FLT 5690 
FLT 570 C 
FLT 5 7 1C 
FLT 572 C 
FLT 5730 
FLT 574 C 
FLT 575C 
FLT 576 C 
FLT 5770 
FLT 57SC 
FLT 5790 
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TREA*FTREE. READDS 


1 - 


SUBROUTINE R EADDS ( NSR t T LD » TL tR 0 . S U PRO 6 * T L CS * TLS tIRBt I SI 

READ 

1C 

2 

C 


READ 

20 

3 

C 

THIS SUBROUTINE READS THE DORMANT PARAMETERS (FAILURE RATES GR 

READ 

30 

4 

c 

FACTOR) AND THE SWITCHING OPTIONS AND PARAMETERS IF THERE 

READ 

40 

5 

c 

ARE ANY STANDBY BLOCKS IN THE RBD 

READ 

5 C 

6 

c 


READ 

60 

7 


DIMENSION SWFRCBI5C) » TL DS (50 ) « TL S ( 5G ) 

READ 

7 C 

3 


DIMENSION ROt 50 * IS (50). TLD (50) *TL( 5G) • IRB ( 50*:30) 

READ 

30 

S 

c 


READ 

90 

10 

c 

READ DORMANCY PARAMETERS. 

READ 

ICO 

11 

c 

DFACT r DORMANCY FACTOR TO BE MULTIPLIED BY STANDBY ACTIVE LAMBDASREAD 

lie 

12 

c 

IF IDUMMY NON ZERO AND DFACT - C - READ LAMBDA DORMANT - TL0(N8). 

READ 

120 

13 

c 


READ 

13 C 

14 


READ 4 10 1 » IDUMMY. DFACT 

READ 

140 

15 

4101 

FORM AT (12. El 2.7 .F1C. 7 .5 5X *11 ) 

READ 

15 G 

16 


IF( IDUMMY. EQ.O) GO TO 4160 

READ 

16C 

17 


IF (DFACT. GT.L. ) 30 TO 4160 

READ 

17C 

13 

4150 

READ 4101. No. TLAM6D. DUMMY. LAST 

READ 

180 

19 


TLD ( NB)=TLAMED 

READ 

19C 

2C 


IF( LAST.LT.6) GO TO 415C 

READ 

200 

21 


GO TC 4190 

READ 

210 

22 

4160 

DO 4130 J-i.NSR 

READ 

220 

23 


NS-ISl J) 

READ 

23G 

24 


KEN D- IRB ( NS .1 ) +1 

READ 

240 

25 


DO 4170 K-2.KEND 

READ 

250 

26 


N3- IRS( NS. K ) 

READ 

260 

27 


TLD ( NB > = DFACT*TU NE ) 

READ 

27 C 

23 

4170 

CONTINUE 

READ 

2 80 

29 

4180 

C ONTINUE 

READ 

25 0 

3C 

C 


READ 

3 00 

31 

c 


READ 

31 C 

32 

c 

READ SWITCHING OPTIONS AND DATA FOR EACH SENSE BLOCK. 

READ 

320 

33 

c 


READ 

33C 

34 

c 

0 = PERFECT 3 WIT CHlNG ( PROBABILITY SWITCH WORKS EQUALS 1.0) 

READ 

3 40 

35 

c 

1 = CONSTANT PROBABILITY THAT SWITCH WORKS. 

READ 

35C 

36 

c 

2 - OORMANT FAILURE RATE FOR SWITCH. 

READ 

360 

37 

c 

3 r DORMANT AND ACTIVE FAILURE RATE FOR SWITCH. 

READ 

37C 

38 

c 

4 - NOT AVAILABLE 

READ 

3 30 

3 S 

c 

5 r NOT AVAILABLE 

READ 

39 C 

40 

c 


READ 

4 00 

41 

4190 

DC 4 7 CC J-l t NSR 

READ 

410 

42 


READ 4101. Nd . TLAMBD. SPR03. ISOPT 

READ 

4 20 

43 


IF (ISCPT.ES. L) GO TO 46 CO 

READ 

430 

44 


GO TO ( 461G.4 62G.463C. 4640.4650) . ISOPT 

READ 

4 4C 

45 

46CC 

SWPRCB(N6)=1.C 

READ 

450 

46 


TLD S ( N8) = C.C 

READ 

460 

47 


TLSt NB ) -C .0 

READ 

470 

48 


GO TC 4680 

READ 

4 80 

49 

4610 

SWPRC6(NE)-SPRC6 

READ 

490 

50 


TLD S( NB ) -0 • 0 

READ 

500 

51 


TLSt NB ) -C .C 

READ- 

510 

52 


CO TO 4S8C 

READ 

5 20 

--53 

4 62 C- 

SWP_R0b(NB)=1.0 

READ 

530 

54 


TLOS( N3 )-TLAM30 

READ- 

"5 40 - 

55 


TLSt N8)=C.O 

READ 

550 

56 


GO TO 4680 

READ 

5 60 

57 

4630 

SWPROB(Nfc)-l.C 

RE AD 

570 

53 


TLD S( N8)-TLAM3D 

READ 

5 80 

59 


READ 4101. NB.TLAMBD 

READ 

59 C 

60 


TLS(NB)=TLAM60 

READ 

6 00 

61 


GO TO 4 6 8 C 

READ 

610 

62 

4640 

CONTINUE 

READ 

620 

63 


GO TO 468C 

READ 

630 

64 

4850 

CONTINUE 

READ 

6 40 

65 

4680 

TLD ( NB ) "0 .0 

READ 

650 

66 


TL( NB ) -0 • C 

READ 

660 

67 


RO(NBl=l.C 

READ 

670 

68 

4700 

CONTINUE 

READ 

6 80 

69 


RETURN 

READ 

69C 

70 


END 

RE AO 

70 C 
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TREE*FTREE. MODIFY 


1 


SUBROUTINE MO DIFY ( NON A CT , MATC SW, A • IPRINT * NRSW 1 

MODY 

1C 

2 


COMMON /ALLSUfc/MLtM2tM 3 *K4 t M5 1 K 6 *M7 tL 1 * L3 • L4 t L5 t L 6 

MOCY 

2 C 

3 


C0MM0N/SU3/IPATH(2C#5O), IAND( 49* 12) t IANDZ ( 2 5, 1 2 ) • I A N, 3L3, M 8 

MODY 

3C 

4 


DIMENSION NON ACT ( 5C t 20 > . M ATC SW C 5C > • N AC V 1 5C 1 

MOCY 

4 C 

5 


INTEGER GATE, A(3>,3L3 

MODY 

50 

6 


DATA N D/C 110 CC00505D5/KN B CHN/QQC77770G3CGC/ NS / 03CCQ0CC5C 5 C5 / 

MOCY 

6 C 

7 


L9=M1~1 

MODY 

7C 

8 


L8=M8-1 

MOCY 

8 C 

9 


DO ICC 1=1,50 

MODY 

90 

1C 


ICC N AC V (11 = C 

MODY 

ICC 

11 


NRSW = 0 

MODY 

110 

12 

C 


MODY 

12 C 

13 

C 

CONDENSE PATHS BT ELIMINATING SERIES ANOS AND ORStSAVlNG LAST* 

MODY 

13C 

14 

c 


MOCY 

14 C 

15 


4 1C DC 30 IK = 1 , L 9 

MODY 

150 

16 


IF (IPATH(1,IK>-3L3>31,3C*31 

MOCY 

18 0 

17 


31 00 35 1= 1 * L8 

MODY 

170 

18 


32 IF (IPATHtI+1 ,IK ) -BL 3 ) 36 ,30,36 

MOCY 

iec 

19 


36 IF( GATE! IP ATH II, IK) ).NE. GATE< IPAT Hi 1 + 1 # IK) ) ) GO TO 35 

MODY 

19C 

20 

c 


MOCY 

20 C 

21 

c 

CHECK TO SEE IF I 

MODY 

2 1C 

22 

c 


MOCY 

22 C 

23 


37 IF C GATE! IP ATH < I , IK ) >-A < 1 )) 34, 35 , 34 

MODY 

2 3C 

24 


34 DO 33 J=I »L8 

MOCY 

24 C 

25 


33 IPA TH( J,IK)=IPATH(J+1,IK ) 

MODY 

2 50 

26 


IPATH(M8,1K1=6L3 

MOCY 

26C 

27 


GO TO 32 

MODY 

2 70 

28 


35 CONTINUE 

MCCY 

28 C 

29 


3D CONTINUE 

MODY 

2 90 

30 

c 


MOCY 

3DC 

31 

c 

MODIFY GATES WITH DORMANCY. GOOD OR G INPUTS TO STANDBY GATES 

MODY 

3 1C 

32 

c 

ARE CHANGED TO D GATES, AND SWITCH OR S GATES APPENDED TO EVERY 

MOCY 

32 C 

33 

c 

D GATE. THIS ENSURES THAT ALL GATES HAVE AT LEAST 2 INPUTS IN 

MODY 

33C 

34 

c 

ADDITION TO CHANGING THE FORM TO ONE WHICH WILL DIRECTLY CONVERT 

MOCY 

34 C 

35 

c 

TO A RELIABILITY BLOCK DIAGRAM. D GATES ARE TREATED AS AND GATES 

MODY 

3 50 

36 

c 

AND S GATES ARE TREATED AS OR GATES* WITH THE GATES* D GATES 

MCCY 

36 C 

37 

c 

CONNECT IN PARALLEL AND S GATES CONNECT IN SERIES. 

MODY 

3 70 

38 

c 


MOCY 

38 C 

39 

c 

MATCSW IS ARRAY CONTAINING SWITCHES. 

MODY 

3 90 

4C 

c 

NONACT IS ARRAY WITH DORMANT BLOCKS AND RESPECTIVE SWITCHES. 

MOCY 

40C 

41 

c 


MODY 

4 1C 

42 


K 5 = C 

MOCY 

42 C 

43 


DO 528 IK = i * Ml 

MODY 

4 30 

44 


DO 522 I = 1 ,2 G 

MOCY 

44 C 

45 


JGT = IPATHd »IK) 

MODY 

4 50 

46 


JG = GATEIJGT) 

MOCY 

46 C 

47 

c 

IS GATE A S GATE • • * 

MODY 

4 7C 

48 


IF 1 JG.EQ.At6) ) GC TO 512 

MOCY 

4£C 

49 

c 

IS GATE A 3 GATE • # • 

MODY 

4 9C 

50 


IF 1 JG *EQ * A (6 d GO TO 513 

MOCY 

5CC 

51 

c 

IS GATE NOT A 0 GATE * • • 

MODY 

510 

52 


IF C JG.NE.AC7J) 30 TC 522 

MOCY 

52 C 

53 

c 

SET KS .NE. 0 , SINCE THIS IS 0 GATE. 

MODY 

5 30 

54 


KS = 1 

MOCY 

54 C 

55 

c 

REMOVE LETTER PART OF THE GATE. 

MODY 

550 

56 


523 JS = AND (JGT #KNBCHN * 

MODY 

56 C 

57 

CHANGE TO S GATE . 

MODY 

570 

58 


JS = ORUS.NS) 

MODY 

58 C 

59 

C 

TO 8E C OR NOT e£ 0 IS THE QUESTION. 

MODY 

5 90 

60 


IF <KS> 529,527*529 

MOCY 

6CC 

61 


529 DO 560 J - 1. Ml 

MODY 

610 

62 

C 

THIS IS A C GATE# SO FIND THAT S GATE TO WHICH THIS D GATE BELONGS. 

MODY 

62 C 

63 


IF (JS.EQ. IPATHd. J) > GO TO 561 

-MODY- 

-6-3G- 

64 

— 

56G-C0NTTNUE 

MODY 

64 0 

65 


561 NACV(J ) = NACV(J ) + 1 

MODY 

5 50 

66 


MK r N A C V ( J l 

MOCY 

66 C 

67 

c 

PLACE THIS SLOCK WITH A 0 GATE INTO NONACT-WITH ITS RESPECTIVE S GATMODY 

67C 

68 


NGNACT 1 J #MK ) = IK 

MOCY 

68C 

69 


KS = 0 

MODY 

6 90 

70 


GC TC 527 

MOCY 

70 C 

71 

c 

REMOVE LETTER PART OF THE GATE. 

MODY 

7 1C 

72 


513 JG = AND (JGT.KNbCHNI 

MOCY 

72 C 

73 

CHANGE FROM G GATE TC A D GATE. 

MODY 

7 3C 

74 


IPATHd, IK1 = OR ( JG.ND > 

MOCY 

74 C 

75 


GO TO 523 

MODY 

7 5C 
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78 


512 NRSW z NR SU * 1 


MGCY 

76 C 

.77 

C 

this is a_ s gate# so list it in matcsw 

AND ADD IT TO THE NUM3ER 

OF SMODY 

7 7C 

78 

C 

FOUND IN THIS TREE THUS F AR . 


MOCY 

780 

79 


MAT CS W t NftS W ) = IK 


M 00 Y 

7 9C 

80 


GO TO 522 


MCCY 

8CC 

81 

C 

SHIFT ALL GATES FOLLOWING THE D GATES 

ONE POSITION IN THE ARRAY 

AHEMODY 

310 

82 

c 

SO THAT A S GATE CAN BE INSERTED AFTER 

THE D GATE, 

MCCY 

82 C 

83 


527 L = 19 - I 


MODY 

8 3C 

84 


DO £25 N z 1 *L 


MOCY 

SAC 

85 


IPATH(M8 - N ♦ 1 1 IK 1 z I PAT HC H8 - N* 

IK) 

MODY 

8 50 

68 


525 CONTINUE 


MOCY 

sec 

87 


IPATH{I+1*IK) - JS 


MODY 

8 7C 

88 


GO TO 528 


MOCY 

££ C 

39 


522 CONTINUE 


MGDY 

8 9C 

90 


528 CONTINUE 


MODY 

50 C 

91 

c 

OUTPUT NONACT* MATCSW AND THE MODIFIED 

TREE AS A CHECK. 

MODY 

910 

92 


IF INRSW.EG.O > 60 TO 578 


MOCY 

92 C 

93 


IF l IPRINT • EQ» 0 ) GO TO 57 8 


MODY 

3 3C 

94 


DO 578 J z 1 t NR S W 


MOCY 

94 C 

95 


NON z MATCSW< J ) 


MODY 

950 

96 


PRINT 516 » NON 


MOCY 

96 C 

97 


516 FORMAT ( * C S WI T CH • 13 r ' CONTROLS THE FOLLOWING BLOCKS*) 

MODY 

970 

98 


DO 530 K zit 2 C 


MCCY 

98 C 

99 


IF* NONACT! NON f K) .EQ.OI GO TO 531 


MODY 

3 90 

IOC 


53C CONTINUE 


MC C Y ICC C 

101 


531 IZK-1 


M0DY101C 

102 


PRINT 52C » l NON ACT t NON . K) *K z 1,1) 


M0CY1C2C 

103 


520 FORMAT (1H .2CI3) 


MCDY103C 

104 


576 CONTINUE 


MO CY 104 C 

105 

c 



M00Y1C5C 

1C6 


576 RETURN 


M0CY1CS C 

107 


END 


M0DY1C7C 


JPL. Technical Memorandum 33-512 



trea*ftree.prteq 


1 


SUBROUTINE PRTEQ( ITEMPtIPRINT • NlfLL6) 



PRES 

10 

2 

C 




PR EG 

20 

3 

c 

THIS SUBROUTINE SETS UP THE COMPLIMENTS 

OF 

THE EQUIVALENT BLOCKS 

PRES 

30 

4 

c 

IN THE SECOND PART OF THE ITEMP ARRAY. 

IT 

ALSO PRINTS THE 

PR EG 

4 0 

5 

c 

THE EQUIVALENT BLOCKS. 



PREQ 

50 

6 

c 




PR EG 

6 C 

7 


DIMENSION ITE MP ( 41 • 23 ) 



PREQ 

7C 

8 


M6-LLE-1 



PREQ 

80 

9 


IF( IPRINT.EQ.G) GO TO 2041 



PREQ 

90 

1C 

14 0 

PRINT 9G5 



PREQ 

ICC 

11 

905 

FORMAT C1HC17HEQUIVAL ENT BLOCKS) 



PREQ 

110 

12 

2041 

DO 123 I-itNl 



PREG 

120 

13 


IF( IPRINT.EQ.G) GO TO 2342 



PREQ 

130 

14 


PRINT 90S t (I TEMPI J » I ) » J -2 tLL 6 > 



PREG 

140 

15 

306 

FORMAT! •0*40131 



PREG 

15C 

16 

2042 

DC 125 K=2tLL6 



PREG 

160 

17 

125 

ITEMP! K^M6 t I)--irEMP(K»I 1 



PREQ 

17C 

18 

123 

CONTINUE 



PREQ 

18 C 

19 


RETURN 



PREQ 

190 

2C 


END 



PREQ 

20C 


24 
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TREA *FT REE « 

PRR8D 




1 



subroutine: p rredi ie ? mi ) 

PRED 

1C 

2 

C 


_ - - _ — ■ — 

“PR90 

" 20 

3 

C 


THIS SUBROUTINE PRINTS THE RELIABILITY BLOCK DIAGRAM. 

PRED 

30 

4 

C 



PRBD 

40 

5 



DIMENSION 13(15.2*50 

PR BD 

5 G 

6 



PRINT 632 

PR3D 

6C 

7 


632 

FORMAT ( *0 RELIABILITY BLOCK DIAGRAM*) 

PRBD 

70 

8 



DC 18 1 = 1 f Ml 

PRBD 

80 

9 



IF (lElltltl ) )21i22t21 

PRED 

9 C 

ID 


22 

IF( ISUt2t 11)21.18. 21 

PR3D 

ICO 

XI 


21 

K=IB(1f1fI)+1 

PRBD 

11C 

12 



IFCK.LE.X) GC TO 40 

PRBD 

120 

13 



PRINT 202.1. (I3( J.l*I).w-2.K ) 

PRED 

13 C 

14 


2C2 

FORMAT! 1H0 5HB LOCK 13 f 3Xt 3 HINPUT14 ( I3»1X) ) 

PRBD 

140 

15 



GO TO 41 

PRBD 

150 

16 


40 

PRINT 205 • I 

PRBD 

160 

17 


2D 5 

F CRM AT ( 1H D5H BLC CKI3* 3Xi 5H INPUT ) 

PRBD 

17 C 



41 

K=IB(1.2tI»+l 

PRBD 

180 

XS 



IFCK.LE.1) GC TC 42 

PRBD 

19 C 

20 



PRINT 203f (IB(J# 2fI)i J=2tK> 

PR3D 

2 CO 

21 


203 

FORK AT ! 11 XfoMCUTPUT * 14( 13 tlX )) 

PRED 

21C 

22 



GC TO 13 

PRBD 

220 

23 


42 

PRINT 2C6 

PRBD 

230 

24 


206 

FORMAT! 1 IX f6H OUTPUT) 

PRBD 

240 

25 


18 

CONTINUE 

PRBD 

25 0 

26 



RETURN 

PRBD 

260 

27 



END 

PRBD 

270 
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TREA*FTREL.IESTE Y 


1 


SUBROUTINE 16 STd Y f $ • N SR* IS . IR B . 18 # IPRINT *L 3 *L 5 ) 



IBSY 

10 

2 

C 





IB SY 

20 

3 

C 

STANDBY 3L0CKS MUST FOLLOW ALL OTHER INPUT BLOCKS 

IN 

INPUT LISTS 

IBSY 

30 

4 

C 

FOR IB. 




IBSY 

4 C 

5 

C 

SUBROUTINE IBST3Y REARRANGES ANY ELEMENTS IN 18 THAT 

NEED IT* AND 

IBSY 

50 

6 

c 

ALSO PRINTS THE SENSE BLOCKS. 




IBSY 

60 

7 

c 





IBSY 

70 

8 


DIKE NSIQN NHCLO( 15 > * IS ( 5C > * I RB ( 50 * 30 1 * 1 6 ( 15 * 2 * 5C 1 



IBSY 

8 C 

9 

c 





IBSY 

90 

1C 


00 2C2C I-ltNSR 




IBSY 

ICC 

11 


INS-IS ( I ) 




IBSY 

110 

12 


CO 2200 J-l • L3 




IBSY 

12 0 

13 


NHOLO ( J ) =0 




IBSY 

130 

14 

2200 

CONTINUE 




IBSY 

14 0 

15 


NUM-C 




IBSY 

150 

16 


KEND=IR3CINS»l)*i 




IBSY 

16 0 

17 

c 





IBSY 

170 

18 


DO 2250 L~2 t L3 




IBSY 

18 0 

19 


IF< I3(L*1. INS J.EO.O) GO TO 226C 




IBSY 

190 

2C 


DO 2240 K-2 . KEND 




IBSY 

200 

21 


IF< I3<L*l*IN51.NC.IR31lNS*K)) GO TO 224C 



IBSY 

210 

22 


NUMr NUM + 1 




IBSY 

22 C 

23 


NHOLD ( N U M ) - IR B < IN S # K ) 




IBSY 

2 30 

24 


IBtL *1.INSI=C 




IBSY 

240 

25 


GO TO 2250 




IBSY 

2 50 

28 

2240 

CONTINUE 




IBSY 

260 

27 

2250 

CONTINUE 




IBSY 

270 

28 

2260 

IF(NUM.NE.C) 60 TO 220C 




IBSY 

280 

29 


PRINT 22 70* XNSiNUMflNS 




IBSY 

2 90 

3C 

2270 

F CRM AT ( * 1 ERR OR IN IBST6Y SUB. AT 2270. 

*/*DERROR DURING SENSE BLOCIBSY 

30 0 

31 


1 K f 1 3 1 * INPUT REARRANGE TO STORE ACTIVE 

BLOCKS BEFORE 

STANDBY 3L0CKIBSY 

310 

32 

2 S IN 1 6 ( * 1 3 t ' 1 1 • * 13 * ’ 




IBSY 

320 

33 


RETURN 1 




IBSY 

3 30 

34 

C 





IBSY 

34 0 

35 

C 

SHIFT X6< ftd»i ' INS) TO LEFT ELIMINATING 

ZEROES. 



IBSY 

3 50 

36 

c 





IBSY 

360 

37 

2300 

ILEFT-C 




IBSY 

370 

38 


00 235C J-2 t L3 




IBSY 

380 

39 


IF( IS C J»1*INS ).NE.C> GO TO 23 30 




I8SY 

3 90 

4G 


I LEF T-ILEFT+1 




IBSY 

400 

41 


GO TO 2356 




IBSY 

410 

42 

233C 

IF(ILEFT.EO.O) GO TO 2350 




IBSY 

420 

43 

C 





IBSY 

4 30 

44 

C 

MOVE NB TC LEFT. 




IBSY 

440 

45 

C 





IBSY 

4 50 

46 


Mr J-ILEFT 




IBSY 

460 

47 


I B ( M* It INS 3U.lt INS) 




IBSY 

4 70 

48 

2350 

CONTINUE 




IBSY 

48 0 

49 

C 





IBSY 

4 90 

50 

C 

ADD NB IN NHOLD AT END. 




IBSY 

50 0 

5 1 

c 





IBSY 

510 

52 


KbEG=L3-lLEFT*l 




IBSY 

520 

53 


Mr C 




IBSY 

530 

54 


DO 237C KrKB EG • L3 




IBSY 

540 

55 


M=M«-1 




IBSY 

5 50 

56 


IB(K*l*INS)rNHOLD(MJ 




IBSY 

560 

57 

2370 

CONTINUE 




IBSY 

5 70 

58 

C 





IBSY 

58 C 

59 

C 

PRINT SENSE BLOCK AND STANDBY BLOCKS OF 

THAT SENSE 

BLOCK. 

IBSY 

5 90 

6C 

C 





IBSY 

60 C 

61 


IF( IPRINT.EQ. C) GO TO 2020 




IBSY 

610 

62 


DC 2C25 Kr2 t L5 




IE SY 

62C 

S3 


IF( IR0( INS *K) . EQ.D) 30 TO 2026 




IBSY 

6 3 G_ 

6j4 — 

65 

, — - -2 0 2 5 
2026 

-CONTINUE— ' 

Kr K-l 




IB SY 
IBSY 

64 C 
6 50 

66 


FRINT 2021 » INS * < IRE ( INS t J ) » jr2,K) 




IBSY 

66C 

87 

2021 

FORMAT { * OS ENS E SWITCH * 13 * • CONTROLS BLOCKS * * 2 «/ 1H 

2914) ) 

IBSY 

6 70 

68 

2020 

CONTINUE 




IBSY 

sec 

69 


RETURN 




IBSY 

6 90 

70 


END 




IBSY 

70 0 


26 
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TREA*FTREE.SUE7 


1 



SUBROUTINE SUS71IB) 

SU37 

10 

2 

C 


- __ — - - - - ------ - - - - ■ - 

- SUE7 

- 20 

3 



DIMENSION NSAV£C43)tI3Cl5f2f5G) 

SU67 

30 

4 



CCMMON/ALLSUfc/Ml «M2 *M3 * K4 *M5 *M6*H7*LI*L3*L4*L5*L6 

SU 67 

4 C 

5 



COMMON/ SUS/ IP ATH(20t5G)t IAND(49t 12* • IANDZ* 25* 12 ) 1 1 A N* BL 3t M 8 

SU37 

50 

6 



INTEGER BL3 

SUE7 

80 

7 



L9 : Ml - 1 

SU37 

70 

8 



DC 13 I = 1 * LS 

SUE7 

80 

9 


13 

NSAVE( I)=C 

S US 7 

90 

1C 

C 



SUE7 

ICC 

11 



N-0 

SUB 7 

110 

12 

c 



SUE7 

12 C 

13 

c 


THE 142 LOOP TAKES THE SLOCK ENTRIES STORED IN I AND THAT 00 NOT 

SU87 

130 

14 

c 


HAVE THEIR INPUTS SPECIFIED AND PUTS THESE BLOCKS IN NSAVE. 

SUE7 

14 C 

15 

c 



SUB7 

150 

16 


7 

DO 142 K - 2 1 L 9 

SUE7 

16 C 

17 



IF( IANDIK. I AN ))179. 142*179 

SUB 7 

17C 

18 


179 

L-IANO(KtlAN) 

SUB7 

18 0 

19 



IF ( IB (2.1* L) >142*180*142 

SUB7 

19C 

EC 


180 

N = N4l 

SUE7 

2CC 

21 



NSAVE 1 N)zL 

SUB 7 

210 

22 


142 

CONTINUE 

SUE7 

22 0 

23 

c 



SU37 

2 30 

24 

c 


GET RID OF DUPLICATE ENTRIES IN NSAVE 

SU £7 

24 C 

25 

c 


THE 157 LOOP IS NEEDED BECAUSE THERE IS NO GUARANTEE THAT I AND 

SUB 7 

250 

26 

c 


DOES NOT CONTAIN DUPLICATE ENTRIES* 

SUE7 

26C 

27 

c 



SU87 

2 70 

28 



N2=N-I 

SU E7 

280 

29 



DO 157 K=I.N2 

SUB 7 

2 90 

30 


85 

IF (NSAVE (K) ) 148 *147 *148 

SUE7 

30 C 

31 


143 

Kl-K+i 

SUB 7 

310 

32 



DC 158 K2-K1 * N 

SUE7 

32 0 

33 



IF(NSAVE(K2)-NSAVE(KiJ158*159*l53 

SUB 7 

330 

34 


158 

CONTINUE 

SUE7 

34 0 

35 



GO TO 157 

SUB 7 

350 

36 


15 S 

DO 84 K 2 - K * N 2 

SUE7 

36 C 

37 


84 

NSAV£(K2)-NSAVE(K2+1) 

SU37 

370 

38 



NS A VE ( N ) -C 

SUE7 

3 80 

39 



GO TO 85 

SUE7 

39 C 

40 


157 

CONTINUE 

SUS 7 

400 

41 

c 



SU E7 

410 

42 

c 


CONNECT THE OUPUTS OF THE SUBDIAGRAM STORED IN IANDZ TO THE INPUTSSUB7 

420 

43 

c 


OF THE SUBDIAGRAM STORED IN IAND (BUT WHOSE LOOSE ENDS ARE 

SUE7 

43 C 

44 

c 


STORED IN NSAVE). 

SUB7 

4 4C 

45 

c 



SUE7 

45 C 

46 


147 

DO 794 K-2 » L3 

SUB 7 

4 60 

47 



IF (IANDZ (K*I AN) ) 181*7 94 *181 

SUE7 

470 

48 


181 

M-IAND21 K* IAN ) 

SUB 7 

4 8C 

49 



DC 160 K2=ltM8 

SUB7 

49 0 

50 


160 

I P A TH ( K2*M)-6L3 

SUB7 

5GC 

51 



DO 143 K1 -1 • N 

SUE7 

510 

52 



LM- NS AVE t K1 > 

SUB7 

5 2G 

53 



IB(K1*1»2 *M)=lM 

SUE7 

530 

54 


143~ 

IBCKniLMriM 

- SUB 7- 

5 40 

55 


794 

CONTINUE 

SUE7 

550 

56 



RETURN 

SUB7 

560 

57 



END 

SUE7 

570 
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TREA*FTREE.GATE 


1 

C 




GATE 

ID 

2 

C 

♦ ♦♦EXPLANATION CF FORTRAN CONTROL CARD*** 



GATE 

2 C 

3 

c 

A COMPILER ERROR IN THE UNIVAC 1103 FORTRAN LEVEL 

7 E COMPILER 

GATE 

30 

4 

c 

INHIBITS RETURNING THE CORRECT VALUE IN GATE. 

WHILE N 

WILL 

GATE 

4 G 

5 

c 

HOLD THE CORRECT VALUE* THE STATEMENT G ATE" N HILL NOT 

TRANSFER GATE 

50 

6 

c 

THIS VALUE TO GATE. THIS HAS 5 EEN FIXED ON THE 

LEVEL 

8 A 

GATE 

ec 

7 

c 

COMPILER* THUS THE CONTROL CARD. NOTE THAT THE 

8 A WILL NOT 

GATE 

70 

8 

c 

BECOME THE SYSTEM STANDARD BECAUSE OF INC CMF ATI 

BILITY 

WITH 

GATE 

8 C 

9 

c 

ENCODE/ DECODE. IF MAKING MODIFICATIONS TO THIS 

PROGRAM# DO 

NOTGATE 

90 

1C 

c 

USE ENCCDE/DECGDE WITHOUT REWRITING THIS SUBROUTINE TO 

WORK 

GATE 

ICC 

11 

c 

WITH THE LEVEL 7E COMPILER. 



GATE 

110 

12 

c 




GATE 

12 C 

13 


INTEGER FUNCTION OATE(IP) 



GATE 

130 

14 


COMMON/GT/BLl 



GATE 

14 0 

15 


INTEGER 3L1 



GATE 

15D 

1C 


IJ-IP 



GATE 

ICC 

17 


N-BLi 



GATE 

17C 

18 


FLD(C*C*N)=FLD(C*£tlJ) 



GATE 

180 

19 


GATE=N 



GATE 

190 

20 


RETURN 



GATE 

20 0 

21 


END 



GATE 

210 


28 
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TREA*FTREE .SETIB 


- 1 - - 



SUBROUTINE _ SET-13 l-IB-tMItL 31- - 


-- -- — ----- 

SET 

- 10 

2 

C 




SET 

20 

3 

C 

THIS SUBROUTINE PUTS THE QUANTITIES 

OF 

INPUTS ANO OUPUTS TO EACH 

SET 

30 

4 

c 

BLOCK INTO THE FIRST ELEMENT OF THE 

IB 

ARRAY. 

SET 

40 

5 

c 




SET 

50 

6 


DIME NSICN IE U5»2»50 1 



SET 

6C 

7 


DO 1021 J - 1 • Ml 



SET 

70 

8 


IF { (I8(2«i»J> ♦ 13 ( 2 • 2 • J I 1 • E G . C I 

GO 

TO 1021 

SET 

80 

9 


DC 1G 10 I - 2 f L3 



SET 

90 

ID 


IF (IB(IfltJ) .EG. 0) GC TC 1C1I 



SET 

ICO 

11 

1C1C 

CONTINUE 



SET 

110 

12 

1011 

IB (ltliJ) - 1—2 



SET 

120 

13 


DO 1C 12 I - 2 * L3 



SET 

130 

14 


IF ( 18 ( 1 1 2 i J I . EQ. Cl GC TO 1C13 



SET 

140 

15 

1012 

CONTINUE 



SET 

150 

Id 

1013 

IB(1.2tJI = 1-2 



SET 

160 

17 

1021 

CONTINUE 



SET 

170 

18 


PRINT 2C06 



SET 

180 

19 

20 06 

FORHATCirti) 



SET 

190 

ZC 


RE TURN 



SET 

200 

21 


END 



SET 

210 
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